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MASS TRANSFER: MEMBRANE 

To model effectively membrane transport phe- 
nomena that pertain to dialysis, two general 
types of membrane models can be considered: 
homogeneous and porous. Homogeneous mem- 
branes are thought of as structureless continua, 
which, in the case of dialysis, consist of poly- 
mer-liquid gels. Porous membranes are viewed 
as an impervious polymer phase, interpene- 
trated by liquid-filled pores. 

The various classes of "real" membranes de- 
scribed in the preceding section can be divided 
into the categories of homogeneous and porous, 
although this distinction is not always an 
obvious one. Since the range of pore sizes in 
available membranes is relatively continuous, 
no sharp delineation exists between homoge- 
neous and porous. Bean (1972) has suggested a 



Mass Transfer in Irregular Geome- 
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criterion for making this distinction, based on a 
comparison of the diffusive and hydraulic per- 
meabilities of a given membrane. 

Modeling of transport in homogeneous 
membranes is well matched by a phenomeno- 
logical approach and a solution-diffusion 
theory. This theory assumes convective trans- 
port to be negligibly small. The transport prop- 
erties of porous membranes can be described by 
continuum hydrodynamic theory, a mechanistic 
approach that considers the structure of the 
membrane. Hydrodynamic theory takes both 
diffusive and convective transport into account. 
Linear, nonequilibrium (irreversible) thermo- 
dynamics constitutes a more general phe- 
nomenological approach that can be applied to 
transport in either class of membrane. This 
approach is especially useful when multiple sol- 
utes and/or driving forces are involved. 
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Transport in Homogeneous Membranes: 
Solution-Diffusion Theory 

In general, the driving force for mass transfer is 
a potential gradient across the membrane thick- 
ness. Phenomenological models used to de- 
scribe mass transfer processes in homogeneous 
membranes are concerned with the relationship 
between this driving force and the resulting 
flux. In dialysis, the primary flux is that due to 
diffusion, which is driven by the concentration 
gradient component of a gradient in chemical 
potential. While other driving forces may exist, 
leading to additional fluxes (e.g., pressure gra- 
dient component of chemical potential gradient 
causing bulk flow), dialysis is fundamentally a 
diffusive process and is modeled as such by 
solution-diffusion theory. 

The cornerstone of solution-diffusion theory, 
Pick's first law, is a simple phenomenological 
model that describes the potential/flux relation- 
ship for diffusive mass transfer in binary sys- 
tems. It states that the diffusive flux is pro- 
portional to the gradient in concentration, mole 
fraction, or chemical potential at constant pres- 
sure. In the strictest sense, a system consisting 
of solute, solvent, and membrane is ternary 
rather than binary. However, a pseudobinary 
approach can be taken in which solute A dif- 
fuses through component B, where B is defined 
as both membrane and solvent. The molar flux 
of solute A can be described by the following 
general form of Fick's first law: 



N Av = x A (N Av + N Bv ) - c t D AB ^-, 

dy 



(12-1) 



where N Ay and N By are the molar fluxes of A and 
B in the y direction, x A is the mole fraction of 
solute A, c t is the total concentration of A and B, 
and Dm is the binary diffusion coefficient. The 
first term on the right represents bulk flow 
resulting from unequal counter diffusion of A 
and B t whereas the second term gives the diffu- 
sive flux. 

In reasonably dilute solutions, the bulk flow 
term may be neglected and the total molar con- 
centration c t may be considered constant. Now 
Eq. (12-1) can be simplified to 
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FIGURE 12-1. Typical concentration profile for 
diffusive transport across a membrane. In this example, 
the solute is less soluble in the membrane than in the 
external phases. 



N Ay = -D. 



AB ' 



dc* 
dy 



(12-2) 



Assuming a linear concentration gradient and 
integrating Eq. (12-2) across the membrane, 



N Ay = 



/ ^ 



(12-3) 



where c lm and are concentrations within the 
membrane at each of its interfaces and / is the 
membrane thickness. 

Figure 12-1 shows a typical concentration 
profile for diffusive transport across a mem- 
brane. In this particular example, the solute is 
less soluble in the membrane than in either of 
the external solutions. Further, c Xs and are 
the concentrations of that solute in the external 
phases at the membrane/solution interfaces. 

A partition coefficient K can be defined as 
the ratio between membrane and external solu- 
tion concentrations at equilibrium: 



K 



(12-4) 



(Note: In general c Im /c u c 2m lc 2 s)- Combin- 
ing Eqs. (12-3) and (12—4) and simplifying by 
eliminating subscripts, 



N = KD&c/L 



(12-5) 
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Equation (12-5) is a form of Fick's first law 
that is particularly well suited to membrane 
problems because external concentrations can 
be readily measured, while concentrations 
within the membrane, as in Eq. (12-3), can- 
not. 

Defining a membrane diffusive permeability 

by 



P m = KD/l (12-6) 



and a membrane mass transfer resistance by 

R m = II KD = l/P mt (12-7) 
Eq. (12-5) can be written as 

N = P m bc = kclR m . (12-8) 



Transport in Porous Membranes: 
Continuum Hydrodynamic Theory 

Yasuda and coworkers (Yasuda, Lamaze, and 
Dcenberry 1968; Yasuda et al. 1969; Yasuda, 
Lamaze, and Peterlin 1971) developed a con- 
cept of dialytic membranes as homogeneous 
water-swollen gels in which thermally induced 
movement of segments of (coiled) polymer 
molecules leaves an interstitial volume free for 
solute transport. They concluded that the per- 
meability characteristics of a highly swollen 
system cannot be described by a single coeffi- 
cient. Values of solute and solvent permeabili- 
ties depend on the conditions of measurement, 
particularly the magnitude of convective flux 
relative to diffusive flux (i.e.^ the transmem- 
brane Peclet number). They concluded that the 
following generalizations were valid: 
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Fick's second law can be derived either from 
Fick's first law (Crank 1975) or the equation of 
continuity (Bird, Stewart, and Lightfoot 1960). 
It describes one-dimensional diffusion in a li- 
quid film in which there are no chemical 
changes and no flow taking place: 



dc 
dt 



D d 2 c 

dy 2 ' 



(12-9) 



Consider a liquid film of thickness 21 that con- 
tains a solute at concentration c 0 . At time t = 0, 
the solute concentration at both surfaces of the 
film is suddenly increased to c x . Integration of 
Eq. (12-9), with appropriate boundary con- 
ditions, can provide concentration profiles 
within the film and fluxes as a function of 
time. 

Colton and Lowrie (198 1) illustrate the solu- 
tion to the unsteady-state diffusion problem 
posed above in the form of dimensionless con- 
centration profiles. The dimensionless group, 
Dr// 2 , is the Fourier number. When Dt/l 2 = 1, 
the film is approaching a new steady-state con- 
centration. A characteristic diffusion time, t D = 
/ 2 /D, can be defined that estimates the 
equilibration time for diffusion across a liquid 
film (or a liquid-filled membrane). 



1. The permeability of a solute, whose size is 
small compared to a calculated membrane 
pore size, is proportional to the degree of 
hydration of the membrane. 

2. Solute permeability coefficients decrease 
exponentially with increasing molecular 
size, when the latter is expressed in terms of 
molecular cross-sectional area. 

3. Solute rejection coefficients change marked- 
ly when the solute cross section approaches 
the dimension of the channels. 



A different approach was used by Klein, 
Holland, and Eberle (1979) for a series of dia- 
lytic membranes that included not only the 
homogeneous, swollen gels studied by Yasuda 
and coworkers, but also glassy polymers having 
porous structures. In describing transport 
through either type of structure, the essence of 
the problem is how to describe the size distribu- 
tion of the channels through which the solutes 
must pass, and the interactions between the 
solute and the channel material. The Yasuda 
approach was based on the statistically predict- 
able formation and disappearance of void 
volumes in the hydrogel. The interaction be- 
tween solute and membrane matrix material is 
ignored. The approach used by Klein and 
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FIGURE 12-2. Plot of 
normalized diffusivity 
(DJD) as a function of 
solute molecular weight for 
a typical dialysis membrane 
(Cuprophan® 150 PM) 
(reprinted from Klein 1977 
with permission). 
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coworkers based on earlier single pore channel 
models (Renkin 1954; Beck and Schultz 1972; 
Verniory et al. 1973) assumes that the mem- 
brane has a pore size distribution that can be 
represented by a single, hydrodynamically 
equivalent pore dimension. Since neither the 
pore fractional cross section A p nor the actual 
pore length l p can be measured independently, 
an experimental method is used to derive their 
ratio. From the Hagen-Poiseuille equation, the 
hydraulic permeability L p is related to the ratio 
A p ll p and the hydrodynamic pore radius r p by 



L p = 4t A 



(12-10) 



where 17 is the viscosity of the solution. The 
value of L p is determined by measuring the 
volume flux of water across the membrane in 
response to an applied pressure gradient. Using 
a diffusional model (Beck and Schultz 1972; 
Vemiory et al. 1973), the ratio A p ll p can be 
related to the diffusive permeability coefficient 
^by 

^WNtf..^ (12 . n) 



a:, 



where q is the ratio of solute radius to pore 
radius, D is the solute diffusion coefficient in 
the solvent, AT, is a power series in q, and D m is 
the effective solute diffusivity in the membrane. 



100 1000 
Test Solute Molecular Weight 



10 000 



For a solute such as tritiated water, whose 
radius is small compared to the pore radius, q is 
approximately equal to zero, K\ = 1, and the 
above equations can be solved simultaneously 
to yield 



4 



(I2r-12) 



Once the single representative membrane 
"pore" size has been determined in this way, the 
permeability of any other solute of known di- 
mensions can be calculated for the same mem- 
brane. In the absence of any specific solute/ 
membrane interactions, such as ionic or 
hydrophobic bonding, this model is useful for 
predicting membrane permeabilities when q < 
0.6. 

On the basis of either model, dialytic trans- 
port clearly decreases as solute size increases 
because of two effects: Increasing molecular 
size implies lower solute diffusivity and in- 
creasing solute sizes produce more interfering 
collisions with the pore walls. The effect is seen 
in Figure 12-2 where the ratio of diffusivity in 
the membrane D m to the solution diffusivity D 
(defined as normalized diffusivity) is plotted as 
a function of solute molecular weight. At the 
low molecular weight end of the plot, the im- 
pedance of the membrane reduces diffusivity to 
50 to 80% of the free solution diffusivity. Mem- 
brane impedance can be defined quantitatively 
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as (1 - normalized diffusivity). At higher 
molecular weights — where solute sizes ap- 
proach the dimensions of the membrane chan- 
nels — the impedance is more than 99%. 

Irreversible Thermodynamics 

The simultaneous presence of both concentra- 
tion and pressure driving forces results in a 
mixed diffiisive-convective transport process. 
When both diffusive and convective fluxes are 
of significant magnitudes, both must be 
accounted for, as well as their interactions. 
Irreversible thermodynamics provides a 
framework for analyzing these complex in- 
teractions between fluxes and forces. Irrevers- 
ible thermodynamics has been applied to mem- 
branes by Kedem and Katchalsky (1958) and 
Katchalsky and Kedem (1962) based on the 
fundamental work of Onsager (1931a, 1931b). 
In addition to the reviews by Katchalsky and 
Curran (1965) and Bean (1972), this topic has 
been discussed from a membrane perspective 
by Bungay (1986) and Spriggs and Li (1976). 

The hydraulic permeability of a membrane, 
L py is expressed in the form of a coefficient of 
proportionality between volumetric flux (super- 
ficial velocity) J v and transmembrane pressure 
difference Ap as shown in Eq. (12-13): 

J v = L p Lp (Air, = 0). (12-13) 

Equation (12-13) applies in the absence of an 
osmotic pressure difference (A tt s = 0). If this is 
not the case, then osmotic pressure differences 
across the membrane will affect the flux. This 
situation is described by Eq. (12-14): 



(12-14) 



where L pD is the coefficient of osmotic flow. A 
second phenomenological equation defines the 
differential rate of transport of solute and sol- 
vent, J D : 



Jd = L Dp kp + L d Att s , 



(12-15) 



where L Dp is defined as the ultrafiltration coeffi- 
cient and L D is the diffusional mobility per unit 
osmotic pressure difference. 



By Onsager's law, L pD = L Dp . Thus, three 
measurements are required to characterize a 
given membrane-solvent-solute system. One 
experiment is suggested by Eq. (12-13): 
measurement of hydraulic permeability L p in 
the absence of an osmotic pressure difference. 

In the case of a selective membrane with 
solute present, a difference in solute and solvent 
fluxes is observed. This exchange flow is de- 
scribed by Eq. (12-16), and it provides a means 
of measuring L Dp : 

Jd = L Dp bp (Att, = 0). (12-16) 

A third experiment can be done with no 
hydrostatic pressure difference, but with solute 
concentration differing across the membrane 
(A 7r 5 > 0). In this case, both an exchange flow 
and an osmotic flow take place. The exchange 
flow defines L D and is described by Eq. (12- 
17): 

J a = L D A-rr s (A/? = 0). (12-17) 

A useful transformation of Eqs. (12-14) and 
(12-15) is given below. Staverman (1951) de- 
fined a reflection coefficient <r: 



<r = -(L pD /L p ). 



(12-18) 



At zero volume flow, the solute flux J s can be 
related to the osmotic pressure by a solute per- 
meability coefficient <o: 

0) = Jj/Att, (J v = 0), (12-19a) 

where a> is related to the membrane permeabil- 
ity P m through the van't Hoff equation. Thus, 



P m = coRT. 



(12-19b) 



With the coefficients defined by Eqs. (12-18), 
(12-19a), and (12-19b), the transformed phe- 
nomenological equations aire (Kedem and 
Katchalsky 1958; Katchalsky and Curran 
1965): 



(12-20) 
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J s = c,(l - o)J v + G>A7I). 



(12-21) MASS TRANSFER: OVERALL 



Equations (12-20) and (12-21) are known as 
the Kedem-Katchalsky equations. 

For systems involving multiple solutes, cou- 
pling occurs in the sense that the chemical 
potential gradient of each solute affects the flux 
of all solutes. Fick's laws, being limited to 
binary systems, do not account for coupling. 
The nonequilibrium thermodynamics approach 
considers that all fluxes depend on all forces, 
thus providing an appropriate theoretical 
framewoik within which to account for cou- 
pling. 

Nonideal Effects 

The plot of reduced diffusivities as a function of 
molecular weight can be misleading when ap- 
plications are projected from such data. On the 
basis of the reduced diffusivity for the examples 
shown, almost no transfer of solutes greater 
than about 5000 dalton would be expected. 
However, in "real- world" devices mass transfer 
occurs in the presence of differential transmem- 
brane pressures. The need to flow feed and 
dialysate solutions past the membrane surfaces 
creates obligatory pressure drops in the respec- 
tive channels. When these are unequal, transfer 
occurs by convective modes as well as by dia- 
lytic transfer. Moreover, the presence of 
nonpermeating species may create osmotic 
pressures that oppose the applied hydraulic 
pressure differences. Thus, the computation of 
mass transfer in a device is often quite a. bit 
more complicated than in a laboratory test cell 
in which flow and pressure conditions can be 
controlled more closely. 

But even under the more or less ideal cir- 
cumstances of laboratory test cells, nonideal 
effects stem from interactions between the feed 
solution and the membrane. For membranes 
made from hydrophobic polymers, the best de- 
scribed phenomenon is solute adsorption within 
the pores. Strong hydrophobic bonding between 
large solutes, such as peptides and proteins, can 
alter the effective pore dimensions of the di- 
alysis membranes. This is reflected in reduc- 
tions of both the hydraulic permeability and the 
dialytic transfer of large solutes. 



In an overall sense, mass transfer in the dialysis 
process involves three phases: feed, membrane, 
and dialysate. In addition to the resistance to 
solute diffusion encountered in the membrane, 
dif fusional resistances in each of the two liquid 
phases result in concentration gradients adja- 
cent to the membrane. Within these boundary 
layers, solute concentration differs from that of 
the bulk solution. The magnitude of the mass 
transfer resistance across such a boundary layer 
is dependent on the thickness of the layer, 
which depends, in turn, on fluid flow parame- 
ters. Mass transfer coefficients are a convenient 
means of quantifying mass transfer rates in the 
case of the two moving fluid phases, feed and 
dialysate, for both laminar and turbulent flow. 
As shown in Eq. (12-22), the mass transfer 
coefficient k is defined by an expression of the 
same form as Fick's first law [see Eqs. (12-5) 
and (12-8)]. The concentration difference Ac is 
that between the membrane surface and the bulk 
liquid: 



N = *Ac, 



(12-22) 



where k replaces the permeability P m , which in 
Fick's first law is based on a constant film 
thickness. Mass transfer coefficient k is ex- 
pressed in the same units as P mt that is, linear 
velocity, and, like P mf is the reciprocal of mass 
transfer resistance, but unlike P m , is based on a 
variable film thickness. The combined effects 
of simultaneous convective and diffusive mass 
transfer that occur in permeable ducts are diffi- 
cult to analyze theoretically. Mass transfer 
coefficients are useful, both as a means of ex- 
pressing analytical results and as a measurable 
experimental parameter in cases for which anal- 
ysis is inadequate. 

Consider a dialyzer of unspecified geometry 
operating countercurrently. The feedside flow 
is assumed to be laminar, while the dialysate 
flow is assumed to be turbulent. The velocity 
profiles for this case are illustrated by Keller 
(1973), as shown in Figure 12-3(a). Solutes in 
the feed are convected parallel to the membrane 
surface, diffuse through the membrane phase, 
and are then transferred by convection again in 
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FIGURE 1Z-3. (a) Velocity profiles for the case of 
laminar feed flow and turbulent dialysate flow, (b) 
Concentration profiles corresponding to the velocity 
profiles in (a) (after Keller 1973). 



the dialysate. Figure 12-3(b) shows the result- 
ing (qualitative) concentration profiles at a sin- 
gle point with respect to the length of the di- 
alyzer. 

The overall mass transfer resistance R Q can 
be expressed as the sum of the resistances 
associated with each phase: 



R 0 = R f + R m + R d , (12-23) 



where the subscripts /, m, and d refer to feed, 
membrane, and dialysate, respectively. 

The overall mass transfer coefficient k a is the 
reciprocal of R a : 



Vk 0 = 1/^+ \IP m + \lk d . (12-25) 



Considering now the length coordinate, both 
momentum and concentration boundary layers 
develop along the length of a flow conduit as 
shown in Figure 12-4. Most practical devices 
operate such that the momentum boundary layer 
is fully developed within a negligibly short 
axial distance, while, in the same devices, con- 
centration boundary layers develop throughout 
their length. As a result, the local mass transfer 
coefficient, as defined above, varies throughout 
the length of the device and is difficult to mea- 
sure experimentally. A length-averaged, overall 
mass transfer coefficient is a readily measured 
parameter. It expresses the relationship between 
the average molar flux N and a length-averaged 
concentration difference Ac: 



N = kjic. 



(12-26) 



Evaluating k Q can be done analytically for lami- 
nar flow in tubes and rectangular ducts. For 
irregular geometries and/or turbulent flows, 
empirical correlations can be used. In either 
case, the heat transfer literature serves as a 
basis, and dimensionless groups provide an 
efficient means of expressing the results. 

Mass Transfer in Permeable Ducts with 
Laminar Flow 

In some dialyzer designs, in particular, hemodi- 
alyzers for which extracorporeal blood volume 
must be minimized, the feedside flow is chan- 
neled through relatively thin ducts and is lami- 
nar in nature. As stated above, this situation 
lends itself to an analytical approach. 

The Sherwood number (analogous to the 
Nusselt number in heat transfer) is a form of 
dimensionless mass transfer coefficient. It is 
defined by 



Sh = kL/D, 



(12-27) 



k 0 = \IR 0 



(12-24) 



where L is a length characteristic of the duct 
geometry and D is the diffusion coefficient of 
the solute species under consideration. 
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FIGURE 12-4. Flow conduit within a diffusive mass transfer device showing 
development of momentum and concentration boundary layers (after Colton and Lowrie 
1981). 



A length-averaged, feedside Sherwood num- 
ber Sbf can be defined as 



Sly = kfL/D. 



(12r-28) 



This Sbf can be related to experimentally 
measurable concentrations as well as estimated 
analytically for laminar flow in regularly 
shaped ducts. Equations (12-29) and (12-30) 
define a dimensionless concentration c* and a 
dimensionless axial length z* for laminar flow 
in a porous tube: 

c* = (c b - c f )l{c d - Cy), (12-29) 

z* = zDlvd%, (12-30) 

where c b is the average bulk solute concentra- 
tion in the feed duct, cy and Cd are solute con- 
centrations in the feed (as it enters) and di- 
alysate, respectively, z is distance in the axial 
direction, D is the solute diffusion coefficient, v 
is average axial velocity, and d ti is the diameter 
of the tube. 

Equation (12-31) defines a wall Sherwood 
number (Sh w ), the ratio of the concentration 
gradient in the fluid adjacent to the wall to that 
in the wall, and Eq. (12-32) defines a length- 
averaged overall Sherwood number (ShJ in 
terms of the previously defined feedside and 
wall Sherwood numbers: 



1/Sh* = 1/Shn, + 1/Sh/. (12^32) 

Solution of the convection-diffusion equation 
with appropriate boundary conditions gives val- 
ues for Sh 0 , from which values for Slymay be 
calculated using Eq. (12-32). Colton and 
coworkers (Colton and Lowrie 1981; Colton 
1969; Colton et al; 1971) have reviewed various 
solutions for tubes and iectangular ducts. The 
classical Graetz problem is stated for lami- 
nar flow heat transfer in a tube with constant 
temperature at the wall, analogous to con- 
stant wall concentration in the mass transfer 
problem. L6veque's solution (L£veque 1928) 
assumes a linear velocity profile near the 
wall and is valid in the concentration entrance 
region. This solution is a useful approxi- 
mation in many cases and is given in Eq. 
(12-33): 



Sh r = 1.62z*-° 



33 



(12-33) 



Sh w = k^AfJD, 



(12-31) 



Colton (1969) solved the case of a constant 
flux boundary condition with finite wall per- 
meability for rectangular channels, which is 
more pertinent to a typical dialysis process. 
Davis and Parkinson (1970) solved the equiv- 
alent case for tubes. These solutions are shown 
graphically by Colton and Lowrie (1981) and 
Colton (1988). Their graph of the Davis and 
Parkinson solution for tubes is reproduced in 
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FIGURE 12-5. Davis and 
Parkinson solution for laminar-flow 
mass transfer in a tube with 
permeable walls— constant flux 
boundary condition. The family of 
curves represents values of the wall 
Sherwood number Sh* (after 
Colton and Lowrie 1981). 



Figure 12-5, which shows a family of curves 
covering a range of wall Sherwood numbers 
from Shn, = 0.1 to 20. In the concentration 
entrance region (dimensionless length z* < 
10~ 2 ), the L6veque solution is a good 
approximation to the Davis and Parkinson solu- 
tion. At higher values of dimensionless length, 
corresponding to more fully developed bound- 
ary layers, the L6veque solution tends to un- 
derpredict the Sherwood number and does not 
predict the asymptotic value of = 4 predicted 
by both Davis and Parkinson for tubes and Col- 
ton for rectangular channels. 

Subsequently, Cooney, Davis, and Kim 
(1974) and Jagannathan and Shettigar (1977) 
included a variable dialysate-side resistance in 
reformulating the analysis as a conjugated 
boundary value problem for parallel-plate 
and hollow-fiber configurations, respectively. 
These latter approaches require solution by nu- 
merical methods. The analysis of Jagannathan 
and Shettigar for hollow fibers, which also con- 
siders the effects of ultrafiltration, predicts that, 
for constant membrane area, fiber diameter and 
length have very little effect on overall mass 
transfer for low-permeability membranes. 
However, for high-permeability membranes, 
overall mass transfer decreases with increasing 
membrane diameter. 



Mass Transfer in Irregular Geometries 
with or without Turbulent Flow 

The dialysate-side flow in many dialyzer de- 
signs occurs within an irregular geometry and 
may be turbulent as well. Analysis of these 
situations is difficult and the resulting pre- 
dictions tend to exhibit larger errors than in the 
above case of laminar flow in well-defined 
ducts. Instead, correlations from empirical heat 
and mass transfer studies are usually employed. 

One dialysate-side flow configuration that 
has been treated analytically is the shell side of 
hollow-fiber dialysis modules. The starting 
point for these analyses is modeling of laminar 
flow relative to arrays of cylinders. (Most 
hollow-fiber dialyzers contain sufficiently small 
diameter fibers at high enough packing densi- 
ties that shellside flow is laminar). The "free 
surface" or "equivalent annulus" model of Hap- 
pel (1959) applies to flow either parallel or 
perpendicular to the cylinder's axes and is lim- 
ited to low packing densities. Sparrow and 
Loeffler (1959) obtained an analytical solution 
for the parallel flow case that applies at high as 
well as at low packing densities, where it agrees 
with Happel's model. 

Hermans (1978), Noda and Gryte (1979), 
and Gostoli and Gatta (1980) have developed 
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models for hollow-fiber dialyzer performance 
that incorporate dialy sate- side mass transfer re- 
sistance. Hermans makes use of an equivalent 
annulus approximation to model a hollow-fiber 
dialyzer with shellside flow perpendicular to the 
fibers. Gostoli and Gatta apply the same 
approximation to parallel flow in both cocurrent 
and countercurrent directions. They define a 
parameter, a, which is the ratio of the inside 
diameter of the fictitious annulus to the fiber 
outside diameter. Fiber packing density 0 can 
be defined by 



(12-34) 



where N is the number of fibers, d to is the fiber 
outside diameter, and d s is the inside diameter 
of the module shell. The quantity a can also be 
defined as the square root of <I>. For values of a 
< 0.635 (4> < 0.4), estimation of Sherwood 
numbers by the equivalent annulus approxima- 
tion is accurate to within 5% (Gostoli and Gatta 
1980). 

Noda and Gryte (1979) applied the complete 
Sparrow-Loeffler solution to dialysate-side 
mass transfer without the usual simplifying 
assumptions. Their model takes into account 
the interactions between neighboring fibers that 
occur at higher packing densities. Figure 12-6 
shows the values predicted for the various Sher- 
wood numbers by their model. Maxima for both 
dialysate-side Sherwood number, Sh,/, and 
overall Sherwood number, Slv,, are predicted to 
occur at a packing density value of 0.63. Analy- 
ses based on the equivalent annulus approxima- 
tion (Hermans 1978; Gostoli and Gatta 1980) 
neglect interactions between fibers and predict 
Sho to increase with packing density without 
reaching a maximum. 

All of the above authors assume regular 
arrays of cylinders, either triangular or square, 
in their analyses. The authors of the hollow- 
fiber models (Hermans 1978; Noda and Gryte 
1979; Gostoli and Gatta 1980) all recognize that 
real dialyzers represent randomly distributed 
arrays whose behavior can deviate significantly 
from that of ideally distributed models. One 
form of deviation that can occur is channeling 
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FIGURE 12-6. Sparrow-Loeffler solution for 
axial-flow shellside mass transfer, including interactions 
between fibers. Assuming a regular, triangular array, the 
model predicts maxima in both dialysate-side and overall 
Sherwood numbers (Shj and ShJ at a packing density 
of <S> = 0.63 (after Noda and Gryte 1979). 



or maldistribution of flow, which is discussed 
further in Chapter 13. In general, these de- 
viations compromise the efficiency of real di- 
alyzers and the model predictions should be 
viewed as an upper limit of achievable perfor- 
mance. 

Devices that deviate significantly from mod- 
el predictions can be treated with correlations 
obtained empirically for similar devices. This 
approach can be used in both laminar and turbu- 
lent flow and for virtually any geometry. In this 
type of correlation, the Sherwood number is 
usually defined in terms of other dimensionless 
groups whose values can be measured ex- 
perimentally. The Reynolds number, Re, char- 
acterizes fluid flow as a ratio of inertial-to- 
viscous forces, while the Schmidt number, Sc, 
expresses the relationship between momentum 
diffusivity and molecular diffusivity: 



Re = pvL/17, 



Sc = i?/pD, 



(12-35) 



(12-36) 



where p is the fluid density, v is a characteristic 
velocity, L is a characteristic length, 17 is the 
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fluid viscosity, and D is the solute diffusion 
coefficient. 

The general expression for local Sherwood 
number is given by 

Sh = /(Re, Sc, duct geometry). (12-37) 

For a length-averaged Sherwood number, the 
dependency on duct geometry drops out and 
correlations of the form of Eq. (12-38) are 
used: 

Sh = a(Re*, ScO, (12-38) 

where a t b, and c are empirically determined 
constants. 

Spriggs and Li (1976), Yang and Cussler 
(1986), Klein, Ward, and Lacey (1987), and 
Prasad and Sirkar (1988) have reviewed mass 
transfer correlations published for heat and 
mass transfer studies for flow external to tubes. 
Values reported in the literature for constant b 
vary from 0.33 to 0.93. The lower end of this 
range probably results from flows that are lami- 
nar in nature, while the upper end may reflect 
effects such as movement of flexible tubes 
(e.g., hollow fibers) in response to the external 
flow field. A value of 0.67 (%) is generally 
accepted for turbulent flow. Better agreement is 
seen for the value of c, which ranged from 0.32 
to 0.38 and is usually fixed at 0.33 (Va). Con- 
stant a varies considerably and appears to de- 
pend on the geometry of the system. 

If a given system is similar in its geometry 
and fluid flow to one for which an established 
correlation exists, then applying the published 
values for the constants a, b, and c will prob- 
ably predict its performance with reasonable 
accuracy. If this is not the case, then a Wilson 
plot may be used to correlate data from the new 
system. In this procedure, first used by Wilson 
(1915) to analyze heat transfer resistances, the 
relationship_ between dialysate mass transfer 
coefficient k d and dialysate flow is assumed to 
follow a relationship defined by 

h = Ag§, (12-39) 

where Q d is the dialysate flow rate and A and b 
are constants specific to the flow geometry and 



o 




FIGURE 12-7. Wilson plot used to conelate 
shellside mass transfer data from devices for which a 
specific correlation has not been established. 



solute under study. Then, substituting in Eq. 
(12-25) gives 

\ik a = \ik f + VP m + 1/AG5- (12-40) 

The quantity k Q is measured at several values of 
Q d and \lk 0 is plotted versus \IQ d as shown in 
Figure 12-7. A value for b is chosen that gives 
the best straight-line fit to the data. The in- 
tercept of this line coiTesponds to infinite di- 
alysate flow rate at which \lk d is zero. Thus, the 
value of the intercept equals {\lk f + VP m ) and 
l/kd at any Q d can be determined by subtracting 
the intercept from the corresponding value of 
llk 0 . 



Dialysis-Specific Mass Transfer Parameters 

In the previous section, mass transfer was dis- 
cussed for the liquid phases adjacent to the 
membrane. For practical applications, this anal- 
ysis must be extended to describe mass transfer 
for the dialyzer as a whole. Neglecting 
convective contributions, we can write the 
following equation for the moles of solute, dQ it 
transferred across an element of membrane of 
length dz and area dA m per unit time (Figure 
: 12-8): 

dQi = k 0 (c f - c d ) dA m9 (12-41) 
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FIGURE 12-8. Diffusive mass 
transfer across an element of a 
membrane of length dz and area dA„ 
(reprinted from Klein, Ward, and 
Lacey 1987 with permission). 
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where c f and c d are the solute concentrations in 
the feed and dialysate, respectively, in the ele- 
ment of length dz. Equation (12-41) can be 
integrated along the length of the membrane to 
give the following equation for overall solute 
transfer rate: 



Qi = k () A„ 



Ac%, - Ac* =0 | 
ln(Ac%,/Ac% 0 ) 



(12-42) 



where Ac* =0 and Ac*=*i are the transmembrane 
concentration differences at either end of the 
device. For countercurrent flow, Ac'JLq = ty — 
c do and Ac* =1 = c/ 0 - c di . This leads to the 
following equation for overall mass transfer in a 
dialyzer with this flow configuration: 



Qi = k 0 A t 



I iCfp - C di ) - {Cfi - c do ) \ 

\ ln[(c fo - c dl )f(c di - c do )] J 



(12-43) 



A similar expression can be derived for cocur- 
rent flow by redefining Ac* =0 and Ac*-i- 

If ultrafiltration is negligible, inlet and outlet 
flows are equal for both the feed and dialysate 
streams. Using this assumption, the following 
overall mass balances may be written for each 
stream: 

Qi = Qf{c fi - C fo ) = Q d (c do - fifl). (12-44) 

The performance of a dialyzer can be described 
in terms of its dialysance D*, which is defined 
as the rate of mass transfer divided by the con- 
centration difference between inlet feed and in- 
let dialysate (Klein, Ward, and Lacey 1987), 
that is, 



D* = Q < 



(12-45a) 



Clearance, defined by Eq. (12-45b), is de- 
rived from renal physiology and is similar to 
dialysance except that dialysate solute concen- 
tration is not considered (Colton and Lowrie 
1981): 



c = Sl. 

C fi 



(12-45b) 



If the solute concentration is zero in the 
dialysate (as in the case of a single-pass di- 
alysate system), then dialysance and clearance 
have the same value. 

Combining Eqs. (12-43), (12-44), and (12- 
45a), the following expression is obtained for 
the dialysance of a countercurrent dialyzer in 
terms of membrane properties and hydrody- 
namics (*, 
Q d Y 



OJ A m ) and operating conditions (£/, 



D * = Q I exp[Mmd - QflQ*)IQi\ - 1 

' \exp[Mm(l - Q f IQ d )IQf\ - Q f fQ di 



(12-46) 

Note that if the feed and dialysate flows are 
essentially equal (Q f = Q d )> then 



D* = Q f 



k A 



Q/ + k 0 A„ 



(12-47) 



The foregoing assumes that the solute is distrib- 
uted freely in the solvent phase. This assump- 
tion may not always be valid; for example, in 
protein solutions, small solutes may bind to the 
protein molecules and exist in equilibrium be- 
tween the free and bound states. In such cir- 
cumstances, modifications of Eq. (12-46) must 
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be used to determine dialysance (Farrell et al. 
1974). 

Equation (12-46) can be used to estimate the 
degree of separation of two solutes by a given 
dialyzer. For any solute, the extraction ratio E 
of solute from the feed stream is given by 

E = Qrfr-(Qr-g*>» (12 _48) 

For negligible ultrafiltration, this reduces to 
E = c fi- c fo (12 _ 49) 

If the inlet concentration of solute in the di- 
alysate is zero (c di = 0), Eq. (12-49) further 
reduces to 

D* 

E = . (12-50) 

Qi 

The extraction ratio represents the fraction of 
maximum solute concentration change that can 
be attained under a given set of operating con- 
ditions. It is analogous to the effectiveness of a 
heat exchanger. The terms of Eq. (12-46) can 
be combined through definition of the following 
dimensionless parameters: 

N t = k 0 AJQ ft (12-51) 

Z = Q f /Q d . (12-52) 

Here N t is defined as the number of transfer 
units and is a measure of the mass transfer size 
of a dialyzer. Extraction ratio E for countercur- 
rent flow can now be expressed as 

£= l-exp[^(l-2)] / 
1 -ZexpHV f (l -Z)] 

Similar expressions can be derived for other 
flow configurations including well-mixed di- 
alysate flow (Michaels 1966). For cocurrent 
flow, 

E= l-exp[-M(l + Z)1 
1 + Z 



For perpendicular flow, 

E = — I [S n {N,)S n {N,Z)), (12-55) 

where 

S n {Y) = l-exp(-y) I — ). 

The term "perpendicular flow" refers to the case 
wherein the feed and dialysate streams flow at 
90° relative to each other. This term was chosen 
instead of "crossflow" because "crossflow" is 
used interchangeably with "tangential flow" to 
describe flow tangential to the membrane sur- 
face. 

The relationships between E t N r , and Z are 
shown graphically for each of the three flow 
configurations in Figures 12-9, 12-10, and 12- 
11. The separation factor a jk in a dialytic pro- 
cess can be considered to be the ratio of the 
fractional masses of two solutes removed from 
their common feed stream, under a given set of 
operating conditions. By applying Equation 
(12-48) to each solute, we can see that a jk is 
equal to the ratio of fractional extractions of the 
two solutes. Thus, Figures 12-9, 12-10, and 
12-1 1 can be used to estimate the relative sepa- 
ration of two solutes from a knowledge of their 
overall mass transfer coefficients, membrane 
area, and feed and dialysate flow rates. 

Contribution of Convection to Mass 
Transfer 

The above derivations are based on the assump- 
tion that convectiye mass transfer is negligible; 
that is, significant ultrafiltration is not occur- 
ring. If significant ultrafiltration does occur, 
mass transfer will be enhanced by the addition 
of a convective component to the diffusive mass 
transfer described above. Ultrafiltration will 
occur from the feed to the dialysate in response 
to a pressure gradient, • which may be either 
applied, in order to concentrate the feed, or 
obligatory, as a consequence of the geometry of 
the device and the desired feed flow rate. Less 
commonly, if the feed is highly concentrated 
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FIGURE 12-9. Extraction ratio £ as a function of number of transfer units N t in a 
countercurrent flow dialyzer. The graphs of Figures 12-9, 12-10, and 12-11 can be 
used to estimate the relative separation of two solutes in each of three flow con- 
figurations. 
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FIGURE 12-10. Extraction ratio £ as a function of number of transfer units N t in a 
cocurrent flow dialyzer. 
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PERPENDICULAR FLOW 
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FIGURE 12-11. Extraction ratio £ as a function of 
number of transfer units N, in a perpendicular flow 
dialyzer (reprinted from Colton and Lowrie 1981 with 
permission). 

and the dialysate dilute, ultrafiltration may 
occur from the dialysate to the feed in response 
to an osmotic gradient. In the first case, ul- 
trafiltration will enhance solute transport while 
in the second it will impede it. 

As discussed earlier, at a steady state the 
overall solute flux J s at a point on the membrane 
can be considered to consist of the sum of a 
diffusive and convective component: 



J s = P m SLc 5 + JM - <r)c s 



(12-56) 



where the first term on the right side, P m Ac S9 
represents diffusive transfer and the second 
term, J v (l - a)c Si convective transfer (Spiegler 
and Kedem 1966). The average solute concen- 
tration within the membrane, c s , can be es- 
timated using Eqs. (12-57), (12-58), and (12- 
59) (Villaroel, Klein, and Holland 1977): 



Ac, 



(12-57) 



where c w is the feedside solute concentration at 
the membrane wall, which can, in turn, be 
estimated using concentration polarization 
theory, of which Eq. (12-58) is an example 
specific to dialysis: 



+ 0 + (1 - <r)fr d fc f 
c f (1 - cr)(l + 0 + tfr ' 



(12-58) 



where £ = [exp(Pe) - I]" 1 (where Pe is typically 
<0.1 in dialytic processes) and ip = [exp(0) - 



1 } _1 . For hollow-fiber membranes , 0 is given by 
6 = 0.7097J - , (12-59) 

where z is the axial distance from the hollow- 
fiber inlet to the point at which c w is being 
evaluated. 

The relative magnitudes of diffusive and 
convective transfer are functions of the mem- 
brane permeability P m and the reflection coeffi- 
cient a-. For species of small molecular size, the 
resistance to diffusive transfer through a mem- 
brane is low, that is, P m is large, and diffusive 
transfer always significantly exceeds con- 
vective transfer at low Peclet numbers. Howev- 
er, as solute molecular size increases, mem- 
brane permeability decreases logarithmically 
(Farrell and Babb 1973), whereas reflection 
coefficients increase at a much lesser rate and 
the relative importance of convection to over? 11 
mass transfer increases. This effect is illustrate 
in the following example. 

Figure 12-12 shows the membrane per- 
meability and reflection coefficient as a func- 
tion of molecular weight for a cellulosic dialysis 
membrane (Wendt et al. 1979). Consider two 
solutes, A and B t having molecular weights of 
200 and 2000 daltons, respectively, both having 
feedside concentrations at a point on the mem- 
brane surface of 0.1 g/cm 3 . Assuming negligi- 
ble dialysate-side concentrations, Ac s = 0. 1 
g/cm 3 (referring to both solutes). Using Eqs. 
(12-57), (12-58), and (12-59), an estimate of 
the value of c s can be made. For both solutes, 
Eq. (12-58) estimates the value of c w /c f to be 
unity (c w as c f ). Then, Eq. (12-57) gives c s - 
0.067 g/cm 3 for both solutes A and B. 

For solute A, P m = 3.6 x 1CT 4 cm/s and 1 - 
a = 0.8. Substituting these values into Eq. 
(12-56), the following result is obtained: 

J 5 = (3.6 x 10- 5 + 0.054y v ) g/s • cm 2 . 

(12-60) 

As the ultrafiltration rate increases from 0 to 0.5 
x cm/s, J s increases from 3.6 x 10~ 5 
g/s-cm 2 to 3.9 x 10~ 5 g/s-cm 2 ; that is, convec- 
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FIGURE 12-12. Membrane permeability and reflection coefficient as a function of 
solute molecular weight for a typical dialysis membrane (Cuprophan® 150 PM) (after 
Klein, Ward, and Lacey 1987). 



tion increases overall solute transfer by 8%. For 
solute B, P m = 0.40 x 10^ cra/s and 1 - cr = 
0.6. Under these conditions, Eq. (12-56) re- 
duces to 

J s = (0.40 x 1<T 5 + 0.04/ v ) g/s • cm 2 . 

(12-61) 

Now, for an increase in ultrafiltration rate from 
0 to 0.5 x KT 4 cm/s, J s increases from 0.40 x 
10" 5 to 0.60 x 10~ 3 g/s-cm 2 , an increase in 
overall mass transfer of 50%. 

The above considerations deal with a point 
on the membrane. Application of these con- 
cepts to the dialyzer as a whole requires that the 
feed and dialysate solute concentrations at ev- 
ery point on the membrane be expressed in 
terms of known concentrations. Such an 
approach has been developed by Schindhelm, 
Farrell, and Stewart (1977), who by assuming a 
zero inlet dialysate concentration and a linear 
dialysate concentration profile, developed the 
following expression for the dialysance of sol- 
utes greater than 300 dalton molecular weight 



in a countercurrent dialyzer in the presence of 
simultaneous diffusion and convection. [Note 
that the equations for X and Y given in this refer- 
ence contain misprints; the correct expressions are 
as given in Eqs. (12r-63) and (12r-64).] 

X 



1 + Y 



(12-62) 



where 



x M, + <kO— > ;jgg_ tl _ exp(Z)] , 

Qf,-Qur z 



(12-63) 



Q/iQdoZ 2 



x 1 1 + gl.u - «p(z)] 



Qfi 



(12-64) 
(12-65) 
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A similar analysis, but one that requires a nu- 
merical methods solution, has been described 
by Jaffrin, Gupta, and Malbranq (1981). 

The contribution of convection to overall 
mass transfer, illustrated previously for a point 
on the membrane, can now be estimated for the 
dialyzer as a whole. Referring to the previous 
example, consider the contribution of convec- 
tion to dialysance for the two solutes A and B in 
a dialyzer containing 10 4 cm 2 of cellulosic 
membrane operating in a countercurrent con- 
figuration with feed and dialysate flow rates of 
8 and 16 cm 3 /s, respectively. For such a mem- 
brane, the overall mass transfer coefficients of 
solutes A and B are on the order of 3.23 x KT 4 
and 0.37 x KT 4 cm/s, respectively, while the 
reflection coefficients are on the order of 0.2 
and 0.4, respectively. Using these data in Eqs. 
(12-62) through (12-65), we can calculate that 
the dialysance of solute A will increase by 12% 
(from 2.47 to 2.78 cm 3 /s) as the ultrafiltration 
rate increases from 0 to 0.5 cm 3 /s, while the 
same increase in ultrafiltration rate will result in 
an increase of 84% (from 0.36 to 0.66 cm 3 /s) in 
the dialysance of solute B, This example further 
demonstrates the important contribution of con- 
vection to overall mass transfer for larger 
molecular weight species. 

NOTATION 
General Notation 

See the General Notation section at the beginn- 
ing of this handbook. 



Special Notation 





membrane area, L 2 


A p 


total pore opening area/membrane area, 




dimensionless 


C 


clearance, l?lt 


d ti 


inside diameter of the tube or hollow fi- 




ber, L 


d to 


fiber outside diameter, L 


D* 


dialysance, lilt 


E 


extraction ratio, dimensionless 


k 


mass transfer coefficient, Lit 


K 


distribution coefficient, dimensionless 


l P 


pore length, L 



L D diffusional mobility per unit osmotic 
pressure difference, Utp 

L Dp ultrafiltration coefficient, Utp 

L p hydraulic permeability, Utp 

L pD coefficient of osmotic flow, Utp 

N molar flux, moVL 2 t 

N( molar flux of species i, mol/L 2 f 

N t number of transfer units 

P m diffusive permeability of solute in di- 
alysis, Lit 

q ratio of solute radius to pore radius, di- 
mensionless 

Q d dialysate flow rate, I? It 

Q f feed flow rate, I? It 

Quf ultrafiltration rate from the feed to the 
dialysate, I? It 

R mass transfer resistance, tIL 

y distance, normal direction, L 

z distance, axial direction, L 

z* distance, axial direction, dimensionless 

Z ratio of feed to dialysate flow rate, Q/IQ* 
dimensionless 

Greek Letters 

j] viscosity, MILt 

a reflection coefficient defined by Eq. (12- 

18), dimensionless 
a) solute permeability coefficient, mol/L 2 tp 

Subscripts 
d dialysate 
o overall 
s solute 
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MEMBRANE SELECTION 
Membrane Materials 

When a dialysis membrane separates two 
aqueous phases, the problem of material choice 
is relatively simple. Clearly, the membrane 
pores must be filled with water during use, 
so wettability is a primary consideration. As a 
consequence of this requirement, the predomi- 
nant materials for dialysis membranes are rel- 
atively hydrophilic polymers. At one end of the 
spectrum are cellulose and poly(ethylene-co- 
vinyl alcohol) (Eval), and on the other end 
poly (methylmethacry late), the latter prewetted 
by the manufacturer. In between fall cellulose 
acetate with a degree of substitution (DS) of 
2.5, poly(acrylonitrile-co-methallylsulfonic ac- 
id), and other acrylonitrile copolymers. 

The use of poly(ethersulfone) membranes 
(PES) for hemodialysis would appear to be a 
contradiction to the wettability requirement, but 
186 
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many dialyzers made from this material contain 
a small amount of poly(vinylpyrrolidone) 
(PVP) as an alloying agent to the PES. The PVP 
provides the requisite wetting, although there 
may still be very small diameter pores that are 
inactive during dialysis because of the exclu- 
sion of water. 

The membrane in widest use for dialysis is 
regenerated cellulose. It can be produced by a 
variety of means, but two processes dominate. 
The oldest and most widely used process is 
based on producing a copper amine complex of 
cellulose in a strongly alkaline solution (i.e., 
the cuprammonium process). The solution con- 
taining 5 to 8% cellulose is coagulated by dis- 
sociating the complex and removing excess 
copper and alkali. The characteristics of mem- 
branes produced by this process are very high 
wet strengths, very thin membranes (as thin as 
0.0005 cm), and a uniform cross section. A 
second process for producing cellulose mem- 
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branes relies on casting a cellulose acetate 
membrane from a high-temperature solvent. 
The resulting cellulose acetate membrane is 
then hydrolyzed to de-acetylate the polymer to 
the native cellulose. Because cellulose acetate 
polymers are generally lower in molecular 
weight than the cellulose used in the cupram- 
monium process, these membranes have lower 
strength and stability to alkali. 

The cellulose membranes are atypical in that 
they contain a high percentage of water, yet do 
not lose their mechanical integrity. As will be 
shown later, this high water content correlates 
with rapid solute diffusivity within the mem- 
brane. 

When the membrane separates two differing 
and immiscible phases, the material selection 
must reflect which phase is to be the pore wet- 
ting phase (see Chapter 41 for further details). 
TTius, if the dialytic process is to remove an 
organic acid from an organic solvent by extrac- 
tion with an aqueous base, a lipophilic mem- 
brane might be chosen. The organic phase 
would wet the pores preferentially to the 
aqueous phase and transfer would occur at the 
aqueous/membrane interface. The organic layer 
in the membrane would become depleted and 
must then be considered as a major contributor 
to the boundary layer. If, instead, a hydrophilic 
membrane were selected, the membrane would 
be wet with the aqueous phase (see Chapter 41 
for alternative wetting conditions) and transfer 
would occur at the organic/membrane interface. 
Now the salt of the organic acid would accumu- 
late in the membrane structure until it is carried 
away by the aqueous phase. The boundary layer 
problem is then reversed. A knowledge of the 
relative membrane diffusivities, i.e., organic 
acid in organic solvent in the membrane vis-&- 
vis organic salt in the aqueous membrane phase, 
would aid in predicting the preferred choice of 
membrane material. 

In either of the above cases, the membrane 
material would have to be completely resistant 
to the solvent effects of either phase. Absence 
of solubility is not a sufficient criteria. Many 
membrane structures are produced as metasta- 
ble phases that very slowly undergo compaction 



of their structure to more stable forms. If the 
contacting solvent accelerates this process, the 
membrane life may be shortened such that prac- 
tical applications are impossible. 

Finally, the various fabrication options avail- 
able with a given material must be considered. 
Most current dialysis devices are produced in 
the form of hollow fibers or tubes. Materials 
formed by sintering are difficult to produce in 
long lengths such as are required by tubes. 
Thus v the preponderance of membranes avail- 
able today is produced by coagulating a melt, or 
solution, or a suspension via some extrusion 
process. Because diffusion paths should be kept 
short, dialytic membranes are generally less 
than 0.0050 cm thick. Such thin walls coupled 
with a need to wind, or bend, the tubes dictates 
that the materials must have relatively high ten- 
sile moduli. Consideration of these factors, 
together with an understanding of the phase- 
separation mechanisms that lead to semiperme- 
able structures, is what constitutes the art and 
science of "membranology." 

Membrane Structures 

Membrane morphology provides important dis- 
tinctions between dialytic membranes. For sim- 
ple dialysis applications, such as equilibration 
of microsolute concentrations in two aqueous 
phases, the choices are many. Most dialytic 
membranes have pore radii of >15 x 10~ 8 cm, 
which is significantly larger than the Stokes 
radius of a 500-dalton molecular weight solute. 
However, when selective separation of solutes 
is desired, other constraints apply. Dialytic 
membranes with pore radii in the 75 x 10~ 8 - 
100 x lO^-cm range are now available. The 
sieving properties of such membranes allow 
them to retain macromolecules while selective- 
ly removing solutes of "middle" molecular 
weight. These "high-efficiency" dialysis mem- 
branes typically exhibit high levels of protein 
adsorption, which makes a priori prediction of 
selectivity problematic. 

Since device clearance, in the absence of 
significant ultrafiltration, depends on short dif- 
fusion paths within the membrane, the thinner 
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the membrane, the more rapid the mass trans- 
fer. So the first choice would be to make a 
membrane as thin as the material will allow, 
while still providing the requisite mechanical 
stability for potting and pumping pressures. In 
addition to being thin, the membrane must pro- 
vide the required morphology, i.e., porosity 
and pore size. These are often conflicting de- 
mands. High porosity implies reduced load- 
bearing elements in the cross section and, thus, 
limited tensile strengths. Fortunately, the more 
hydrophilic materials, such as cellulose and 
copolymers of vinyl alcohol, can be produced 
as thin, highly efficient hydrogel membranes 
reinforced by crystalline regions. The more 
hydrophobic, glassy polymers, such as poly- 
(sulfones), are produced in asymmetric form to 
provide a thin solute-resistant skin supported by 
a more porous substructure. The thick sub- 
structure provides the mechanical support 
needed. 

The ability to produce membranes composed 
of a solute-retentive skin integral with a more 
porous substructure was a major development 
in membrane technology. This asymmetric 
structure was found to permit production of 
membranes with large L p values that were still 
retentive to small molecules. The major benefit 
of this asymmetry, however, accrues to ul- 
trafiltration and reverse osmosis membranes 
since the substructure does not produce a stag- 
nant boundary layer in such convective pro- 
cesses. With dialysis membranes, the thickness 
of the porous substructure can add significantly 
to the boundary layer problems by producing an 
unstirred layer of fluid. Such layers create un- 
wanted transport resistances. 

Asymmetric membrane structures are pro- 
duced by the coagulation process used to con- 
vert a solution of polymer into a membrane 
structure. The phase-separation mechanisms 
that occur as the polymer is precipitated from 
concentrated solutions include liquid-liquid and 
liquid-gel transitions. The skin is thought to be 
caused by a rapid loss of solvent from the poly- 
mer solution film into the coagulating bath. The 
resulting highly concentrated polymer solution 
precipitates with a morphology different from 
the underlying polymer solution since the in- 



itially coagulated surface acts as a barrier to 
transfer of the solvent out of— and nonsolvent 
into — the lower layers of the solution. For the 
production of asymmetric dialysis membranes, 
the control of such phase-separation mech- 
anisms becomes rather critical because of the 
added demand for thin cross sections in dialytic 
applications. 

Table 13-1 lists some commercially avail- 
able, hollow-fiber dialysis membranes with 
their basic performance parameters. Both gel- 
type and asymmetric membrane structures arc 
represented. Unless otherwise noted, UFR (ul- 
trafiltration rate per membrane area) and P m 
(diffusive permeability) were measured using 
saline solutions. 

MEMBRANE MODULES 
Module Types 
Industrial Dialyzers 

Early industrial dialyzers comprised three 
general configurations: the tank type, the plate- 
and-frame type, and the tube type. Tuwiner 
(1962) reviews these designs as well as those of 
early laboratory dialyzers in considerable de- 
tail. Of the three types, the plate-and-frame 
design was most widely used, while the tube 
type (a bundle of tubular membranes) did not 
achieve commercial importance. 

The Cerini dialyzer, an example of the tank 
design, was developed in Italy in 1928 to re- 
claim sodium hydroxide in the rayon industry. 
It consisted of a 3 x 1.5x1 ,2-m tank, contain- 
ing 50 membrane bags totaling 300 m 2 of mem- 
brane area. Thus, the membrane area per unit 
volume of this dialyzer was 55 m~ l . Each of the 
cotton bags was supported by an internal metal 
mesh frame and was individually connected to 
the water supply. Feed solution flowed into the 
tank at the bottom and exited at the top by 
simply overflowing. 

The patent literature abounds with examples 
of plate-and-frame (or filter press) dialyzer de- 
signs (Tuwiner 1962). By providing improved 
membrane support compared to the tank design, 
these dialyzers allowed the use of much thinner 
membranes. The lower diffusive resistance of 
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TABLE 13-1. Performance Parameters of Some Representative Hollow-Fiber 
Dialysis Membranes. 



UFR 



/VVitamin B l2 



Description 


Material 


(mL/h nr mm Hg) ( 1 (T* cm/min) 


Cuprophan 


/Ta Mi line a 

ueiiuiose 


«*.U J.J 


CI (bnka) 






Cuprophan® 


Cellulose 


3.0 4.0 


D4 (Enka) 






Higjiflux® 


Cellulose 


30 12 


RC-HP400 






(Enka) 






Gambro HF 


Cellulose 


3.8 a - 4.8 


Hospal 


PAN-methallyl sulfate 


27 9.5 


Baxter/Toy oba 


Cellulose acetate 


4.6 5.7 


Fresenius 


Polysulfone 


4.6 5.9 


F-6 






Fresenius 


Polysulfone 


33 17 


F-60 






Toray/Hoechst 


PMMA 


16» 10 


•Values measured from whole blood. 



the thin membranes led to a fivefold to tenfold 
increase in productivity per unit membrane area 
relative to the Cerini dialyzer. On a membrane 
area per unit volume basis, plate-and-frame and 
tank-type dialyzers were comparable. 

One of the more recent examples of the 
plate-and-frame design is the Graver dialyzer. 
This design uses a poly(vinylchloride) mem- 
brane support frame. A "repeat unit" in such a 
dialyzer consists of two such frames — one per- 
fused with dialysate and one with feed solu- 
tion — two sheets of membrane, and the neces- 
sary gaskets. The Graver dialyzer is comprised 
of 150 such units. Figure 13-1 illustrates a 
single-repeat unit in a generic plate-and-frame 
dialyzer. 

In countercurrent operation, the lower den- 
sity dialysate flows downward and the higher 
density feed flows upward in adjacent com- 
partments. Thus, the increasing density of the 
dialysate as solute is transferred to it and the 
corresponding decrease in density of the feed 
both serve to facilitate a balanced flow distribu- 
tion among the cells. Hydrostatic pressure in 
the dialysate cells is maintained at a higher level 
than in the feed cells, which ensures that the 
membranes are pressed against the support 
structure of the feed frames. 



The hollow-fiber dialyzer, developed first 
for hemodialysis, has been adapted for in- 
dustrial applications. Enka AG offers industrial 
dialyzers based on their Cuprophan® regener- 
ated cellulose hollow-fiber membrane. Hollow 
fibers of 200- /xm inside diameter and 16-/im 
wall thickness are interwoven with thread and 
grouped into bundles. Three such bundles are 
contained in the annular space between two 
stainless steel tubes and encapsulated with 
polyure thane at each end, which forms 
tubesheets. The inner tube acts as a conduit 
feeding the shell side of the module, the flow 
exiting at one end of the outer tube. Two such 
submodules of 28-cm effective fiber length are 
joined in series, comprising a module of 22.5- 
m 2 membrane area. The membrane area per 
unit volume of this device is approximately 
2200 wT l , which is a fortyfold increase over the 
tank or plate-and-frame designs. 

While the traditional aqueous/aqueous di- 
alysis process is limited in its industrial applica- 
tions, modern hollow-fiber dialyzers with ex- 
panded chemical resistance are finding new 
application in aqueous/organic solvent extrac- 
tion processes. An example of an industrial 
hollow-fiber device of this kind is produced by 
Sepracor Inc. (more details are available in 



190 IV / Dialysis 



Feed Out 



Dialysate In I 



Dialysate Out <J 
Feed In Hy<wj > 




0 







coil di 
sign w 
brane. 
a coil 
operat 
lophar 
flowec 
plastic 
Kii 
lyzer : 



Dialysate Frame 
. and Separator 



Membrane 



Feed Frame 
and Separator 



Membrane 



FIGURE 13-1. Single-repeat unit from a typical plate-and-frame dialyzer (with gaskets omitted). 
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Chapter 41). These modules were developed 
specifically for applications involving aqueous/ 
nonaqueous extraction processes and are fabri- 
cated from materials that are resistant to organic 
solvents. Available hollow-fiber types include 
regenerated cellulose and poly (aery lonitrile), 
both of 200- /xm inside diameter. The fibers are 
encapsulated in epoxy, forming a cartridge that 
has ends provided with radial O-ring seals. The 
cartridge is installed in a housing, together 
comprising a module. Sepracor modules are 
sized by total available fiber packing volume 
and range from 1.5 mL to 12 L in size (Figure 
13-2). A 12-L module packed with typical re- 
generated cellulose dialysis fibers has a mem- 
brane area of 65 m 2 . Based on overall volume 
occupied, its membrane area per unit volume is 
comparable to the Enka unit described above. 

Hemodialyzers 

As described in Chapter 11, Kolff and Berk 
(1944) demonstrated the feasibility of hemo- 
dialysis and applied the process clinically using 
a "rotating drum" artificial kidney. This device 
consisted of a drum wrapped with cellulose 
sausage casing through which blood flowed 
while the drum was rotated in a bath of di- 
alysate solution. 

In 1956 Kolff and coworkers developed the 
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FIGURE 13-2. Hollow-fiber modules for 
diffusion-based membrane separation processes. The 
modules are sized by total fiber packing volume and 
range from lab scale (1.5 mL) to production scale (12 
L). (Courtesy of Sepracor Inc.). 
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coil dialyzer. Like the rotating drum, this de- 
sign was based on a tubular cellophane mem- 
brane. The tubing was flattened and rolled into 
a coil along with a sheet of plastic mesh. In 
operation, blood flowed inside of the cel- 
lophane tubing around the coil, while dialysate 
flowed axially through the space formed by the 
plastic mesh. 

Kiil introduced a plate-and-frame hemodia- 
lyzer in 1960. Like industrial plate-and-frame 



modules, it was operated in the countercurrent 
flow mode. It was a reusable device whose 
membranes were replaced by disassembly and 
reassembly of the unit. Due to the amount of 
labor required to maintain reusable dialyzers, 
disposable devices of both the coil type and 
plate-and-frame type were introduced (Figure 
13-3). The disposable designs tend to be more 
compact than their reusable counterparts and 
are sold presterilized. Hollow-fiber hemodia- 




FIGURE 13-3. Hemodialyzers of the coil type and plate-and-frame type (Courtesy of 
Baxter Healthcare Corp.; Klein, Ward, and Lacey 1987). 
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FIGURE 13-4. Typical hollow-fiber hemodialyzer (Courtesy of Baxter Healthcare 
Corp.; Klein, Ward, and Lacey 1987). 



lyzers were first introduced in 1966 (Lipps et al. 
1967). The structure of a typical unit is shown 
in Figure 13-4. A bundle of hollow fibers is 
contained in a housing and encapsulated at each 
end forming tubesheets. At each end, a gasket 
and endcap form headers to direct blood flow in 
and out of the lumens of the fibers. Adjacent to 
each tubesheet is a circumferential header, 
which directs dialysate flow in and out of the 
shellside space. The device is geometrically 
similar to a shell-and-tube heat exchanger. 

The housing and endcaps of these modules 
are usually fabricated from a transparent 
engineering polymer such as polycarbonate or 
poly(styrene-co-acrylonitrile). The encapsulant 
that forms the tubesheet is typically a 
polyure thane, and a low-Durometer elastomer 
such as silicone rubber is used for the gasket 
between the endcap and tubesheet. 

The hollow-fiber dialyzer has achieved 
widespread acceptance and has largely sup- 
planted other designs. This is due to the ease 
with which this design can be manufactured, 
its compactness, and its reliable performance 
in extracorporeal systems. While originally 
intended as disposable items, health care 
costs have provided incentive to reuse these 
dialyzers, but only with respect to a single 
patient. 



Module Design 

From a design standpoint, the evolution of di- 
alyzer types followed two distinct lines: those 
intended for industrial applications and those 
intended for medical applications. Specific ex- 
amples from both application areas were given 
in the preceding section. 

Many of the design requirements of utmost 
importance to membrane devices intended for 
medical use have no relevance to industrial 
membrane modules (and vice versa). For ex- 
ample, all hemodialyzers are packaged sterile 
and are intended for single use or reuse only 
with the same patient. Examples of presterilized 
industrial membrane devices exist; however, 
they are the exception and no presterilized in- 
dustrial dialyzers are known to the authors. The 
converse is equally true. An industrial dialyzer 
that utilizes an organic solvent as extractant and 
is to be periodically cleaned with sodium hy- 
droxide must be constructed from chemically 
resistant materials, unlike those used for medi- 
cal devices. 

Scale is another obvious difference. Mem- 
brane devices intended for medical use are de- 
signed to perform mass transfer operations at a 
scale dictated by human anatomy and physiolo- 
gy. Industrial dialyzers, on the other hand, are 
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scaled according to process economic factors 
and are typically much larger than their medical 
counterparts. 

Process economics, in addition to affecting 
scale, also tends to favor long-lived industrial 
dialyzers, while, as mentioned above, hemodia- 
lyzers are intended for single use or limited 
reuse. The economics of health care as applied 
to end-stage renal disease (ESRD) is equally 
important in determining hemodialyzer de- 
sign—no "overdesign" is tolerated in this cost- 
conscious and competitive market. 

To illustrate the process leading to an op- 
timized module design, a hollow-fiber module 
will be designed for a hypothetical dialysis 
process. As described in the preceding sections, 
hollow-fiber modules have become prevalent in 
both industrial and medical applications due to 
their high volumetric efficiency and ease of 
manufacture. 

The hypothetical process for which this 
hollow-fiber module is being designed is an 
industrial process in which the product is a 
compound of 20K-dalton molecular weight (M) 
in aqueous solution. It is desired to separate this 
product from a contaminant whose molecular 
weight is 200 dalton in a continuous process. 
The product solute concentration is 5 g/L and 
required productivity is 50 kg/h. The contami- 
nant concentration is 1 g/L and should be re- 
duced to 0.1 g/L or less by the process. An 
aqueous dialysate solution is to be used, which 
is buffered to pH 7.0. 

Fiber Dimensions and Number 
Yang and Cussler (1986) conclude that the key 
advantage of hollow-fiber membrane modules 
is the high ratio of membrane area to unit device 
volume they provide and that mass transfer 
coefficients for hollow-fiber devices are, in 
general, similar to those for other types of 
equipment. This conclusion provides a clue 
regarding choice of fiber dimensions — 
minimization of fiber diameters (inside and out- 
side diameters) maximizes area-to-volume. 
However, practical lower limits exist for fiber 
diameters. The internal diameter is limited by 
the size of any suspended particles that may be 



present. Wall thickness, which together with 
inside diameter determines outside diameter, 
has a lower limit set by fiber strength require- 
ments. 

For given diameters, fiber length has an up- 
per limit that is determined by pressure drop 
and/or mass transfer parameters. In the case of 
high-flux membranes, excessive translumenal 
pressure drop (TLP) causes high ultrafiltration 
rates, the effects of which are discussed in 
Chapter 12. A further consequence of high TLP 
is the increased likelihood of Starling's flow or 
back-filtration, which can cause a further de- 
crease in separation efficiency. If, at the max- 
imum lumenal flow rate allowable within TLP 
limitations, the desired degree of solute transfer 
occurs in less than one pass of the feed stream, 
then fiber length is excessive. At the other ex- 
treme, if too short a device is designed, ex- 
cessively high flow rates are required to achieve 
a given value of k 0 over a given membrane area 
and recirculation may be required to provide 
sufficient residence time. Another limit on 
"shortness" is manufacturing cost. The shorter a 
device is for a given membrane area, the more 
potential membrane area is lost to waste and the 
more expensive the housing becomes. 

While not all of these optimization parame- 
ters can be included in a single expression, the 
ones that relate directly to fluid dynamics and 
mass transfer performance can be dealt with on 
a quantitative basis. Referring to the hypotheti- 
cal process described above, based on through- 
put alone, the minimum feed flow rate required 
to process 50 kg/h is Q f = 10,000 L/h. Assume 
that the feed solution viscosity is 1 .0 cp and that 
the maximum permissible TLP is 75 kPa (11 
psi), relative to the effective fiber length (ex- 
cluding fiber encapsulated in the tubesheets). 
Equation (13-1) gives the number of hollow 
fibers required to process a given feed flow rate 
within TLP limitations: 

N = %7]LQ f IO:i^TTr%\ (13-1) 

where 17 is solution viscosity, L is effective 
fiber length, Q f is feed flow rate, and r ti is fiber 
internal radius. 
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A further constraint on the number of fibers 
required is the membrane area needed to pro- 
vide a sufficiently high rate of transfer of the 
contaminant. By choosing commercially avail- 
able dialysis fibers covering a range of di- 
mensions, practical values of both fiber di- 
mensions and permeability can be assumed. 
Enka AG offers dialysis fibers of 100, 200, and 
300 /ira inside diameter (dry dimensions), des- 
ignated C3, CI, and D4, respectively. Based on 
the experience of the authors (Klein et al. 1976, 
1977), these fibers will be assumed in the wet 
state to have inside diameters of 125, 250, and 



375 fim and wall thicknesses of 22, 22, and 32 
/xm, respectively. 

To a first approximation, the diffusive per- 
meabilities of all three fibers chosen are the 
same for low molecular weight solutes (<200 
dalton) (Klein et al. 1976; 1977). It is assumed 
throughout the following analysis that passage 
of the 20K-dalton product is negligible through 
the membranes chosen. A convenient way of 
estimating required membrane area (for a first 
iteration of the design) is to refer to clearance 
values published for hemodialyzers of a given 
size. Figure 13-5 gives a value of C = 90 
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FIGURE 13-5. Clearance of various dialysis membranes as a function of molecular weight (after Enka AG 1982). 
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FIGURE 13-6. Number of fibers required as a function of fiber dimensions for 
constant translumenal pressure drop and area. (Used in designing a hollow- fiber dialysis 
module for a hypothetical process.) 



mL/min for a 1-m 2 dialyzer using one of the 
hollow fibers chosen above. Dividing Q f by C 
gives an estimate of required membrane area, 
A m = 1850 m 2 . Equation (13-2) gives the num- 
ber of hollow fibers required to provide a given 
membrane area as a function of fiber di- 
mensions: 



N = A m l{2m ti L)~ 



(13-2) 



Figure 13-6 was plotted to illustrate how the 
above constraints of TLP and membrane area 
limit the design space for the fibers chosen. The 
three curves that originate in the lower left cor- 
ner of the graph were generated by Eq. (13-1) 
using the above-stated values for TLP and 
Q f — each curve represents a different fiber 
inside diameter. Values of fiber number and 
length on or above each of the curves result in 
acceptable TLP values that are below 75 kPa. 
The family of curves that originates in the 



upper left comer of the graph were generated by 
Eq. (13-2) with A m = 1850 m 2 . Values of fiber 
number and length on or above these curves 
result in membrane areas that are sufficient for 
mass transfer. Thus, for the 125-/im-i.d. fiber, 
greater than 10 7 fibers are required for most 
fiber lengths. Hie minimum number occurs at 
the intersection point of the two curves: L = 85 
cm and N - 5.6 x 10 6 fibers. For the 250-pun- 
i.d. fiber, the minimum number isN = 1.0 x 
10 6 at a fiber length of L = 240 cm. The 
corresponding values for the 375-ptm-i.d. fiber 
are N = 3.6 x 10 5 and L = 440 cm. At this 
point in the design process, it would be pre- 
mature to select from among the three fiber 
sizes. A detailed evaluation of mass transfer 
performance may help in the decision. 

Flow Configurations 

Feed and dialysate flows can be either cocurrent 
or countercurrent in the case of parallel flows. 
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Other possibilities discussed in Chapter 12 in- 
clude well-mixed dialysate and perpendicular 
flow. Of these possible flow configurations, 
countercurrent and perpendicular flow offer the 
best potential mass transfer performance. 
Cocurrent operations are indicated when 
minimal convective flow is desired and oblig- 
atory pressure drops are encountered. Under 
cocurrent flow configurations, pressure differ- 
ences are easier to balance across the 
membrane along the length of the device. 
Countercurrent flows are used to maximize 
concentration differences between the inlet di- 
alysate and the outlet feed, where the latter is 
most dilute. 

Based on Eqs. (12-53) through (12-55), ex- 
traction ratio E versus the number of transfer 
units N t has been plotted for three of the four 
flow configurations discussed above. Figures 
12-9 through 12-11 of Chapter 12 show these 
plots of the countercurrent, cocurrent, and per- 
pendicular flow cases. Extraction ratio E is per- 
haps the most meaningful measure of dialyzer 



mass transfer performance. In countercurrent 
flow, at equal values of N t and Z, somewhat 
higher values of £ are obtained than in the other 
two flow configurations. This implies that 
countercurrent flow is the most efficient of the 
three configurations. Recall, however, that N t 
= k 0 A m IQ f . Thus, comparison at equal values 
of N t implies that (if A m and Q f are to be held 
constant) k Q is the same in both countercurrent 
and perpendicular flow hollow-fiber modules 
operating at the same Qj. This is most often not 
the case and k 0 is usually higher for the per- 
pendicular flow case. 

Yang and Cussler (1986) compared pub- 
lished heat transfer correlations for flow ex- 
ternal to tubes in both the parallel and per- 
pendicular flow configurations. They further 
conducted oxygen transfer experiments with 
hollow-fiber modules in which external water 
flow in both flow configurations was tested 
under conditions where membrane resistance 
was negligible. Figure 13-7 compares the pub- 
lished correlations as well as the correlations 
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FIGURE 13-7. Mass transfer and heat transfer correlations for flow external to tubes, 
comparing parallel flow and perpendicular flow. 
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FIGURE 13-8. Relative 
contributions of feed, 
membrane, and dialysate 
mass transfer resistances for a 
typical hollow-fiber dialyzer 
operating in countercurrent 
flow (reprinted from Colton 
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mission). 
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resulting from their experiments. In both cases, 
the Sherwood number at a given value of 
Reynolds number is more than an order of mag- 
nitude higher in perpendicular flow compared 
to parallel flow. 

Consider the case in which membrane mass 
transfer resistance is significant, which is more 
relevant to dialysis. Figure 13-8 shows relative 
contributions of feed, membrane, and dialysate 
resistances as a function of solute molecular 
weight for a typical hollow-fiber dialyzer op- 
erating in countercurrent flow. For a solute of 
200 dalton M, R d contributes about 25% of the 
overall mass transfer resistance. Applying the 
results from Yang and Cussler (1986) for per- 
pendicular flow, a dialyzer containing the same 
fiber as in Figure 13-8, but with dialysate in 
perpendicular flow, should show a 25% de- 
crease in R 0 due to Rd becoming negligible. 

By applying the above example to the plots 
of Figures 12-9 and 12-11 of Chapter 12, k Q 
(perpendicular) is equal to 1.33&* (parallel) at 
the same Reynolds number and comparison 
should be made between the perpendicular flow 
case at a value of N t that is 33% higher than that 
used in the countercurrent flow case. Compared 
in this way, for a value of Z = 0.5, per- 
pendicular flow attains E = 0.9 at N t = 5, 
while, for countercurrent flow at the same 
Reynolds number, N t = 3.75 and E = 0.9. 
When enhanced dialysate mass transfer in per- 
pendicular flow is taken into account, per- 



pendicular flow and countercurrent flow are 
predicted to give comparable overall mass 
transfer performance in a typical dialytic pro- 
cess. 

Shellside Pressure Drop 
Another consideration in choosing between per- 
pendicular and countercurrent flow is shellside 
pressure drop. The above mass transfer com- 
parison was done at equal Reynolds numbers, 
which is an appropriate basis for a pressure drop 
comparison as well. Happel's (1959) free- 
surface model predicts pressure drops for lami- 
nar flow relative to arrays of cylinders. He 
states that his model shows good agreement 
with experimental results (on arrays of cylin- 
ders) up to Re = 100. Referring to the Reynolds 
number definition, Eq. (12-35) of Chapter 12, 
the Reynolds number for shellside flow is 
obtained by inserting fiber outside diameter for 
the characteristic length L in Eq. (12-35) of 
Chapter 12. The wet outside diameters for the 
three hollow fibers chosen are 170, 295, and 
440 /Ltm. Thus, at Re = 100, average velocities 
in the shell for the three fiber sizes are v = 59, 
34, and 23 cm/s, respectively. The free-surface 
model should be applicable since actual hollow- 
fiber modules are likely to operate below these 
velocities. 

Application of the free-surface model (Hap- 
pel 1959) to a moderately packed (4> = 0.5) 
triangular array of the 295-/xm-o.d. fibers gives 
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resulting values for Darcy's constant of K D = 
1 . 1 8 x 10 - * cm 2 for perpendicular flow and K D 
= 3.06 x IQr 6 cm 2 for parallel flow in the same 
array [K D = 7/v/(Ap//)]. Since K D is inversely 
proportional to pressure drop per unit length, 
the model predicts a pressure drop that is 2.6 
times as large for perpendicular flow as for 
parallel flow over the same length. However, 
most practical hollow-fiber module con- 
figurations have axial lengths that exceed their 
transverse flow dimension by factors greater 
than 2.6:1. Thus, at this level of analysis, com- 
parable pressure drops would be expected for 
parallel and perpendicular flow. 

Keep in mind that the free-surface model 
applies to rigid, regularly arranged cylinders 
and that real hollow fibers may deviate signifi- 
cantly. Hermans (1978) compares results pre- 
dicted by the free-surface model for per- 
pendicular flow with an empirical study using 
hollow fibers (Dandavati, Doshi, and Gill 
1975) and concludes that measured pressure 
drops exceeded those predicted by the model by 
a factor of 5: 1. Deviations in the direction of 
higher-than-predicted pressure drops can be at- 
tributed to packing of fibers in response to flow, 
deformation of fibers due to the pressure gra- 
dient, nonuniform outside diameter, and "kinki- 
ness." The first two of these causes apply 
specifically to the perpendicular flow case. 
Packing of fibers can be controlled by weaving 
them into a mat or providing local support with 
supplementary fibers. "Kinkiness," which may 
increase pressure drop in parallel flow, can be 
of benefit in perpendicular flow to control pack- 
ing. In general, the pressure drop behavior in 
parallelflow is more predictable than in per- 
pendicular flow, and simple fiber bundle con- 
figurations (e.g., without weaving or supple- 
mentary fibers) give good performance. 

For the reasons stated above, the parallel 
flow regime will be chosen for the module be- 
ing designed. Predictability, simplicity of de- 
sign, and module-to-module consistency are 
good reasons to choose parallel flow for a di- 
alysis process. However, for processes in which 
one of the fluid boundary layers contributes 
most of the overall mass transfer resistance, 
perpendicular flow should be considered. For 



example, for the experiments conducted by 
Yang and Cussler (1986), both membrane re- 
sistance and gas-phase boundary layer resis- 
tance were quite small compared with liquid- 
phase boundary layer resistance. In this case, 
perpendicular flow gave overall mass transfer 
coefficients that exceeded those for parallel 
flow by more than an order of magnitude, as 
described in the preceding section. 

Prediction of Overall Mass Transfer 
Performance 

Use of Figures 12-9 and 12-1 1 of Chapter 12 to 
estimate values for the extraction ratio requires 
knowledge of a dialyzer's overall mass transfer 
coefficient k 0 . An experimentally determined 
value of k Q for the solute, membrane, and di- 
alyzer configuration of interest is to be pre- 
ferred. Estimates of k Q can be made based on 
knowledge of the geometry and fluid dynamics 
of the dialyzer and the nature of the membrane 
and solute. Methods have been described in 
Chapter 12 _and this chapter that allow estima- 
tion of P m , kf, and k d for a variety of situations. 
Regardless of whether an experimental mea- 
surement or a calculated estimate of k 0 is used, 
it is important to make the determination at the 
projected operating temperature because of the 
strong temperature dependence of solute dif- 
fusivity. 

As discussed earlier, the overall mass trans- 
fer resistance (R 0 = \lk Q ) can be considered the 
sum of the individual boundary layer and mem- 
brane resistances; that is, as shown in Eq. (12- 
25) of Chapter 12: 

llk 0 = lik f + \IP m + \lk d . 

If the individual mass transfer_coefficients can 
be estimated, then a value for k Q can be derived 
from Eq. (12-25) of Chapter 12. Since three 
examples of hollow-fiber membranes have been 
selected, published data can be used as a predic- 
tion of P m . Such studies have been done by 
Klein et al. (1976, 1977). By interpolation from 
the results of these studies, P m for 125- and 
250-/im-i.d. fibers, both of which have 22-/xm 
wet wall thicknesses, is estimated to be 4.0 x 
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ltT 4 cm/s for a solute of 200 daltons M and the 
375-/xm-i.d. fiber, which has a wall thickness 
of 32 fim, has an estimated permeability of 3.1 
x KT 4 cm/s. 

"Before attempting to predict k df a value for 
packing density should be chosen. Noda and 
Gryte (1979) predicted $ = 0.63 as the op- 
timum value for packing density, but they, as 
well as Yang and Cussler (1986) and Hermans 
(1978), warn of the effect of channeling on 
shellside mass transfer. If shellside flow is 
poorly distributed among the fibers, then rela- 
tively stagnant regions exist that severely com- 
promise mass transfer. Yang and Cussler 
observed this effect in a close-packed module 
that exhibited flow-independent mass transfer at 
1/10 the rate of the same fibers in a loosely 
packed module. The value of O at which 
channeling becomes a problem is dependent on 
module design factors, such as fiber bundle 
configuration (woven or unsupported), shell in- 
side diameter and ratio of fiber length to outside 
diameter. Of course, it is desirable to approach 
* = 0.63 to maximize both surface-to- volume 
ratio of the module and mass transfer coeffi- 
cient per unit flow rate. The authors' experience 
has been that packing densities in the 50 to 60% 
range can be achieved without undue channel- 
ing provided the fiber length is large in propor- 
tion to both fiber and shell diameters. The ini- 
tially proposed fiber lengths arrived at in Figure 
13-6 certainly qualify. Therefore, <E> = 0.55 
will be chosen as a reasonable value for the 
design at hand._ 

To estimate k d , a correlation developed by 
Wald, Kessler, and Lopez (1990) will be used. 
The basis for this correlation was a series of 
experiments using a number of identical 
hollow-fiber modules with Lid much larger than 
hemodialyzers, which tends to minimize 
channeling. Mass transfer performance of this 
series of modules showed little module-to- 
module variation. Equation (13-3) gives the 
correlation 

Shj = 0.025 Re 0 * 4 Sc 033 . (13-3) 

This correlation is similar to that reported by 
Yang and Cussler (1986) for their loosely 



packed hollow-fiber modules. Equation (13-3) 
predicts higher values for the Sherwood num- 
ber, especially for large values of Lid, than does 
the Yang and Cussler correlation. The value of 
exponent b is similar in both correlations (0.94 
and 0.93) and is difficult to explain in relation 
to laminar flow in bundles of rigid tubes. 
However, this value is typical of shellside mass 
transfer in hollow-fiber modules (Colton 1988). 
The fact that this exponent is much higher than 
predicted or observed for laminar flow in well- 
defined ducts and higher even than that 
observed for turbulent flow can be explained if 
we consider the freedom that hollow fibers have 
to move laterally and orient with the flow field. 
As fluid velocity increases in parallel flow, the 
resulting forces on the fibers tend to improve 
the uniformity of fiber distribution, which, in 
turn, minimizes channeling and improves over- 
all mass transfer.^ 

In estimating k d for the module being de- 
signed, a value of 1100 is estimated for the 
Schmidt number (based on a diffusion coeffi- 
cient of 9 x I0r* cm 2 /s for a 200-dalton solute). 
The Reynolds number will be calculated for 
each of the three "minimum fiber number" 
points selected in Figure 13-6, assuming a feed 
flow rate of 10,000 L/h and a dialysate flow rate 
of 20,000 L/h. Then a value for SYi d is calcu- 
lated using Eq. (13-3) and k d is calculated using 
Eq. (12-27) of Chapter 12. The results of these 
calculations are tabulated in Table 13-2. 

Now the fluid-side mass transfer coefficient 
based on logarithmic mean concentration differ- 
ence can be estimated. First, inserting the above 
estimates of membrane permeability and di- 
alysate-side mass transfer coefficient into Eq. 

TABLE 13-2. Values of Reynolds Number and Di- 
alysate-Side Sherwood Number and 
Mass Transfer Coefficient Corre- 
sponding to Minimum Fiber Numbers 
of Figure 13-6. 



Fiber Outside 




Sh rf 




Diameter (jun) 


Re 


k d (cm/min) 


170 


9.0 


2.0 


0.064 


295 


29 


6.0 


0.11 


440 


55 


11 


0.13 
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TABLE 13-3. Values of Wall Sherwood Number and 
Dimensionless Length Used to De- 
termine Feedside Sherwood Numbers. 





Fiber Inside Diameter (/im) 


Sh w 


z* 


I 


125 


0.40 


1.2 




250 


0.91 


0.61 


A*"" 
y ■ 


375 


1.1 


0.40 



(13-4), values for can be calculated for each 
of the fiber diameters: 



l/k w = \IP m + \lk d 



(13-4) 



Then, an estimate can be made for the wall 
Sherwood number, as defined by Eq. (12-31) 
of Chapter 12, at each fiber inside diameter. 
Values for dimensionless length z* can be 
calculated for the minimum fiber number con- 
ditions for each of the three fibers. The above 
results are tabulated in Table 13-3. 

Referring to Figure 12-5 of Chapter 12, val- 
ues of feedside (lumen-side) Sherwood num- 
bers can be determined from the curves. Feed- 
side mass transfer coefficient kf can then be 
calculated. Overall mass transfer coefficients 
can now be determined by adding each of the 
mass transfer resistances previously determined 
(reciprocal of each k). These values can be 
found in Table 1S-4. 

By comparing the values of membrane per- 
meability and the various mass transfer coeffi- 
cients as resistances (their reciprocals), the rela- 
tive contributions to overall resistance can be 
obtained, as shown in Table 13-5. 

To reiterate, the initial design compares the 
three different fibers at equivalent feed and di- 
alysate flow rates, membrane area, and TLP, 
operating in countercurrent flow. The mem- 



brane area was chosen by approximating overall 
mass transfer using clearance data for a dialyzer 
containing one of the fibers chosen. Feed flow 
rate was chosen as the minimum flow required 
to provide the required throughput of product, 
and dialysate flow rate was arbitrarily set at 
twice the feed flow. From Figure 13-6, con- 
ditions were chosen such that the minimum 
total fiber number for each diameter was 
selected. 

The result, as presented in Table 13-5, is 
favorable in that for all fibers selected, the sum 
of feed and dialysate resistances contributes 
only one-third of overall resistance, the rest 
being contributed by the membrane. 

Further insight into the performance of this 
design and its optimization can be gained by 
referring to Figure 12-9 of Chapter 12, the plot 
of extraction ratio versus the number of transfer 
units for countercurrent flow. Calculating val- 
ues of N t for the three fiber diameters can now 
be done using the k 0 values in Table 13-4. 
These values have been referred to a value of Z 
= 0.5 in Figure 12-9 of Chapter 12 and the 
results are tabulated in Table 13^-6. 

These extraction ratio values predict that the 
initial design conditions will result in between 
67 and 75% removal of low molecular weight 
contaminants, depending on the membrane 
selected. Since the goal is 90% removal, some 
adjustments must be made in module design 
and/or operating conditions to achieve a higher 
extraction ratio. 

At this point, the 375-/im-i.d. fiber will be 
dropped from consideration due to its lower 
level of performance (which is due primarily to 
its thicker wall and resulting lower membrane 
permeability). With the two remaining fibers, 
simply increasing dialysate flow rate will in- 



TABLE 13-4. Feedside Sherwood Numbers and Mass Transfer Coeffi- 
cients and Overall Mass Transfer Coefficients. 



Fiber Inside Diameter (/im) 


Sty 


kf (cm/min) 


k 0 (cm/min) 


125 


.4.3 


0.19 


0.016 


250 


4.3 


0.093 


0.016 


375 


4.3 


0.062 


0.013 
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TABLE 15-5. Contributions of Feed, Membrane, 
and Dialysate Resistances to R 0 in Ini- 
tial Design. 



Hber Inside Diameter 
(/on) 


*/(%) 






125 


8 


67 


25 


250 


17 


68 


15 


375 


21 


69 


10 


crease extraction 


ratio. An 


increase 


from 



20,000 to 40,000 Uh will approximately dou- 
ble the_values of k d associated with each fiber 
(since k d depends on Re 0 * 94 ). Referring again to 
Figure 12-9 of Chapter 12, a further advantage 
results from increasing Q dt reflected in the 
choice of Z = 0.25, which gives higher extrac- 
tion ratio values for a given value of N t . The 
new values are shown in Table 13-7. 

While the changes have brought perfor- 
mance closer to the goal of E = 0.9, further 
increases in Q d will not achieve the goal. In- 
stead, membrane area will be increased by in- 
creasing the total number of fibers, while hold- 
ing fiber length constant. The value of Q f will 
remain unchanged, resulting in a decrease in 
TLP and in kf. If the total fiber number is 
increased by 50%, A m will increase to 2775 m 2 
and TLP will decrease from 75 to 50 kPa. Since 
feedside mass transfer is already operating in 
the asymptotic regionof Figure 12-5 of Chapter 
12, the decrease in kf will be negligible. Di- 
alysate flow rate will be increased by an addi- 
tional 50% to maintain the same Reynolds num- 
ber as in the previous example. Since Q f has not 
increased, Z = 0. 17 in this example. The result- 



TABLE 13-6. Number of Transfer Units N t and Ex- 
traction Ratio E for the Initial Design 
with Z « 0.5. 

Rber Inside Diameter 

0«n) N t E 

125 1.8 0.75 

250 1.8 0.75 

375 1.4 0.67 



TABLE 13-7. Number of Transfer Units N t and Ex- 
traction Ratio E Reflecting a 2x In- 
crease in ft, and Z = 0.25. 



Fiber Inside Diameter 






(>im) 


N t 


£ 


125 


2.0 


0.82 


250 


1.9 


0.81 



ing values of N t and E (from Figure 12-9 of 
Chapter 12) are given in Table 13-8. 

Physical Design of Membrane Modules 
Up to this point, a number of hollow fibers 
operating in parallel have been considered for 
design purposes, without regard to modulariza- 
tion. Since the total number of fibers is in ex- 
cess of 10 6 fibers for the process considered 
above, multiple modules will be required. For 
the two fibers still under consideration, the total 
fiber numbers are now 8.4 x 10 6 and 1 .5 x 10 6 
for the 125- and 250-/xm fibers, respectively. 
Considering the different fiber lengths (85 cm 
for 125 ftm and 240 cm for 250 /tin), total 
module volumes can be calculated as well. At 
4> = 0.55, for the 125-/xm fiber, an internal 
module volume of 295 L is required, while the 
design based on the 250-/Ltm fiber requires 445 
L. Assuming a maximum shell inside diameter 
of 25 cm, seven 85-cm-long modules contain- 
ing the 125-/im fiber or four 240-cm-long mod- 
ules containing the 250-/un fiber would meet 
the requirements. While four of the longer mod- 
ules occupy more total volume than seven of the 
shorter modules, manifolding of four modules 
is simpler. Another practical advantage of the 



TABLE 13-8. Number of Transfer Units AT, and Ex- 
traction Ratio E Reflecting a 1.5x In- 
crease in A m , a 3 x Increase in Q d , and 
Z = 0.17. 



Fiber Inside Diameter 

(*un) N f E 



125 3.0 0.92 

250 2.9 0.91 
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250-/wn-i.d. fiber is that it is more resistant to 
clogging of the lumen. 

Configuration of the modules could be sim- 
ilar to a shell-and-tube heat exchanger or a 
hollow-fiber hemodialyzer (Figure 13-4). 
Shellside flow distribution to a 25-cm-diameter 
bundle of 295-/im-o.d. fibers requires special 
attention, however. Simply providing a baffled 
circumferential header at each end of the shell 
would leave an area of stagnation in the center 
of the bundle. Instead, flow will be introduced 
through a distribution tube at one end and in the 
center of the bundle and will exit circum- 
ferentially at the opposite end. The Enka and 
Sepracor industrial hollow-fiber modules de- 
scribed above use examples of this approach. 

Materials Selection 

For the process described above, the feed is 
relatively noncorrosive (pH 7.0). Material 
selection may be determined more by the 
method used to clean the modules in this case. 
For cellulose membranes, a sanitizing solution 
such as glutaraldehyde or peroxyacetic acid is 
likely to be used. Either 3 16 stainless steel or an 
engineering polymer such as polysulfone can be 
used to fabricate the module housing. The en- 
capsulant used to form tubesheets at either end 
of the module could be a polyurethane, silicone 
rubber, epoxy, or a thermoplastic such as 
polypropylene. 

In the situation where either the feed or di- 
alysate is an organic solvent solution, stainless 
steel is a good choice for the housing, although 
more solvent-resistant polymers such as nylon 
can also be used. The list of possible encapsu- 
lants should be shortened to epoxy or polypro- 
pylene in this case. Elastomeric seals are lim- 
ited to elastomers specifically selected for a 
particular solvent, perfluoro elastomers, or 
fluoro-polymer-coated seals. 

If a membrane is selected that allows steam- 
ing-in-place (SIP), then only stainless steel 
should be considered for a housing material 
with epoxy or, possibly, a thermoplastic as en- 
capsulant. Ethylene-propylene copolymer 
elastomers (EPR or EPDM) are a good choice 
for elastomeric seals in contact with steam. 



Cleaning with caustic solutions requires similar 
materials choices to SIP. If sodium hypochlo- 
rite is to be used as a cleaning agent, then 
stainless steel is not a good choice for a housing 
material. One of the engineering polymers 
would be preferred. 

PROCESS AND SYSTEM DESIGN 
Batch versus Continuous Operation 

Dialysis can be used as a unit operation in two 
basic configurations: as a batch operation or in 
continuous mode. Figure 13-9 depicts a 
schematic of batch dialysis. A reservoir of 
volume V and feed solute concentration c f is 
dialyzed continuously with an inlet concentra- 
tion equal to c di . The feed and dialysate streams 
enter the dialyzer with flows of Q f and Q d , 
respectively. The feed stream exiting the di- 
alyzer returns to the reservoir while the spent 
dialysate is discarded to enter another recovery 
process. Such a configuration could be used for 
stripping low molecular weight contaminants 
from a high molecular weight product. It could 
also be used to recover a low molecular weight 
product by catalytic cleavage of a high molecu- 
lar substrate in the feed reactor. In that event the 
spent dialysate is processed to recover the prod- 




where 
A voli 
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feed 
neglij 
to 



DIALYSATE 

FIGURE 13-9. A batch dialysis process (reprinted 
from Klein, Ward, and Lacey 1987 with permission). 
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uct and not discarded. The former represents, in 
essence, the practice of hemodialysis with the 
reservoir representing the body pool of water. 

A mass balance may be written for the reser- 
voir: 



dt 



(13-5) 



where D* is the dialysance defined previously. 
A volume balance may also be written for the 
reservoir: 



Dialysis reactors may also be operated in a 
continuous mode. This can be represented by 
adding a feed line to the schematic of Figure 
13-9. 

Hie addition of a high molecular weight sub- 
strate at a flow rate Q r to an enzymatic reactor 
whose products permeate the membrane would 
exemplify one case. As the substrate is con- 
verted to dialyzable product at a rate G, the 
increasing concentration must be incorporated 
into the previous equations. Equation (13-5) 
now becomes 



V - V 0 - Q*U 



(13-6) 



(13-11) 



where V 0 is the initial reservoir volume and Q uf 
is the rate of ultrafiltration from the feed to the md me integrated form incorporates the solute 
dialysate. Substituting for V in Eq. (13-5) leads regenerat jon rate and volume expansion as fol- 



to 



(Vb - QufO - QnfC/ - - D *( c / - c <«>- 

dt 

(13-7) 

Rearrangement and integration lead to the 
following relationship for the reservoir concen- 
tration as a function of time: 



lows: 
c/U) 



D*?* + G | t ( (Vo-e^)| '| 



D*-Q Hf + Q r 



(13-12) 



where 




c f (t) = 



+ c 0 | (Vo ~ Qtf» J ^ (13-8) 

where b = (D* - Q^)IQ uf and cj is the initial 
feed solution concentration. If ultrafiltration is 
negligible, that is, = 0, Eq. (13—5) reduces 
to 



dc t 



8 = 



D* ~ Quf + Qr 
Qtf-Qr 



V^- = -D*(Cf- C di ), 03-9) 

dt 

and the reservoir solute concentration is given 
by 



c f {t) = c ?exp -^|+ c di 



1 - exp 



-D*t 



(13-10) 



Staged Operation 

To improve the selectivity of a dialytic separa- 
tion beyond what can be achieved in a single 
stage, a cascade process consisting of several 
stages can be used. Rautenbach (1986) has re- 
viewed the cascade process option as it applies 
to membrane-based processes in general. Noda 
and Giyte (1981) described a multistaged sys- 
tem with two cross-current dialyzer cascades 
and a reverse osmosis unit acting as a nonselec- 
tive solvent stripper. This approach increases 
the slope of the curves in Figures 12-9, 12-10, 
and 12-11 of Chapter 12, relating £ to N t9 
resulting in improved relative separation of two 
solutes of differing molecular weight. Assum- 
ing Qf = Qd throughout the system, the extrac- 
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tion ratio for the system is given by Eq. (13- 
13): 



1 + A7 



(13-13) 



where n is the number of stages in each cas- 
cade. 

To illustrate the efficacy of this approach in 
improving selectivity, assume that two solutes, 
A and B % are being separated in a countercur- 
rent, single-stage dialysis process with Q f = 
Q d . Under these conditions, their extraction 
ratios are£ A = 0.33 and E B = 0.60 with = 
0.5 and N m — 1.5. The resulting selectivity is 
the ratio of E B to E At which is 1.82. In a 
three-stage cascade system, Eq. (13-13) pre- 
dicts E values of 0. 1 1 and 0.77 for a selectivity 
of 7.0. 

NOTATION 
General Notation 

See the General Notation section at the begin- 
ning of this handbook. 

Special Notation 

membrane area, L 2 

initial feed solution concentration, mol/ 
L 3 

clearance, I? It 
dialysance, I? It 
extraction ratio, dimensionless 
conversion rate of the substrate to the 
dialyzable product, mol/r 
mass transfer coefficient, Lit 
Darcy's constant, L 2 
number of hollow fibers 
number of transfer units 
diffusive permeability of solute in di- 
alysis, Lit 

dialysate flow rate, L 3 lt 
feed flow rate, L 3 lt 

ultrafiltration rate from the feed to the 
dialysate, L 3 lt 

inside radius of hollow fiber, L 
mass transfer resistance, tIL 
translumenal pressure . drop, p 



C 
D* 
E 
G 

k 

K D 

N 

N t 

Pm 

Q* 
Qf 
Qu f 

r« 
R 

TLP 



UFR ultrafiltration rate per membrane area, 
Lltp 

initial reservoir volume, L 3 
distance, normal direction, L 
distance, axial direction, L 
distance, axial direction, dimensionless 
ratio of feed to dialysate flow rate, Q f l 
Q d , dimensionless 



V 0 

y 

z 

z* 
Z 



Greek Letters 

ry viscosity, MILt 

O packing density of a hollow-fiber mod- 
ule, dimensionless 

Subscripts 
d dialysate 
o overall 
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General Notation 



Units are given in terms of physical quantities, 
length (L), mass (A/), time (r), temperature (T), 
amount of substance (mol, mole), electric cur- 
rent (A, ampere), electric potential (V, volt), 
energy (£ = ML 2 /t 2 ), and pressure (p = Ml 
Lt 2 ). Special notation is also present in some 
chapters. 



A 
A 



B 

5„ 



d 

df 



activity, various units or dimensionless, 
or particle radius, L 
area, L 2 
species A 

plasticization coefficient for species i, 
L 3 /mol or L 3 /Af 

interaction coefficient between species i 

and j 9 L 3 /mol or I?IM 

channel height, L, or Langmuir affinity 

constant, Lt 2 IM or p~ l , or the ratio of 

frictional force of solute to the bulk 

solution, dimensionless 

species B 

interaction coefficient between species j 
and j, dimensionless 
concentration, mol/L 3 [or cm 3 (STP)/ 
cm 3 polymer] or MIL? 
concentration of species A at the inside 
radius of membrane, mol/L 3 [or cm 3 
(STP)/cm 3 polymer] or Mil? ' 
concentration of species A at the outside 
radius of membrane, mol/L 3 [or cm 3 
(STP)/cm 3 polymer] or MIL 3 
Langmuir sorption capacity, L 3 (STP)/ 
L 3 polymer [cm 3 (STP)/cm 3 polymer] 
heat capacity, L 2 lt 2 T or EIMT 
reflection constant, 1.14 for air, di- 
mensionless 
diameter, L 

diameter of filter fibers, L 



d h hydraulic diameter, L 

d t inside diameter, or molecular diameter 

of species i, L 
d Q outside diameter, L 
d p pore diameter, or particle diameter, L 
d s stirrer or impeller diameter, L 
d t tube diameter, L 
D diffusion- coefficient, L 2 /t 
D dimensionless diffusivity, i\Dlr 2 r wo 
D 0 infinite dilution diffusivity, or charac- 
teristic diffusivity, L 2 lt 
Da Damkohler number, k f d^ l l 2 ID (n = 
order of reaction), dimensionless 
diffusion coefficient for solute-carrier 
complex, L 2 lt 

Brownian diffusion coefficient, L 2 lt 
effective diffusion coefficient, L 2 lt 
diffusivity in the internal phase, L 2 lt 
diffusivity in the membrane phase, L 2 lt 
shear-induced diffusion coefficient , 
L 2 lt 

thermodynamic diffusion coefficient, 
L 2 lt 

energy, ML 2 lt 2 or E, or activation ener- 
gy, ML 2 lt 2 mol or is/mol, or enhance- 
ment factor (F - 1), dimensionless 
activation energy, ML 2 lt 2 mol or £7mol 
Fanning friction factor, or inertia lift 
function, dimensionless 
fugacity of component i, MILt 2 or p 
facilitation factor (ratio of flux with car- 
rier present to flux without carrier), di- 
mensionless, otF = Off/Do, convenient 
dimensionless group in the dual-mode 
sorption model 
F Faraday's constant (9.652 x 10 4 amp • 

s/g-equivalent), A//g-equivalent 
g gravitational acceleration, Lit 2 



Dab 

D b 

O eff 

D< 

D m 

D s 

D T 



/ 

fi 
F 



xvn 



xviii General Notation 



g c gravitational conversion factor, di- 
mensionless [or ft • lb/(lb force * s 2 ) in 
the fps system] 

G Gibbs free energy, ML 2 /t 2 or E 

GPU gas permeation unit, 10"* cm 3 (STP)/ 
(cm 2 • s • cm Hg) 

H enthalpy, ML 2 lt 2 or £ 

H Henry's law constant (c = Hp), t 2 mol/ 
L 2 M [or cm 3 (STP)/(cm 3 polymer • 
atm)], or t 2 IL 2 

H 0 channel half height, L 

&H V latent heat of vaporization (per unit 
mass), L 2 lt 2 or EIM 

i electric current density, A/L 2 

lam limiting current density, A/L 2 

/ electric current, A 

j mass transfer or permeation rate, mol/f 
(or kg-equivalent/s) 

j k mass transfer rate due to diffusional 
transport, mol/r (or kg-equivalent/s) or 
hilt 

j<t> mass transfer rate due to breakage , mol/r 

(or kg-equivalent/s) or Mlt 
J flux, moly L 2 t or MIL 2 t 9 or volumetric 

flux, Lit 

(J) area-averaged volumetric flux, Lit 
Ji flux of species /, moV L 2 t or MIL 2 t 
(J L ) area-averaged volumetric flux from in- 

ertial lift model, Lit 
(j s ) area-averaged volumetric flux from 

shear-induced diffusion model, Lit 
J v volumetric flux, Lit 
k thermal conductivity, MLIt 3 T or E/LtT, 

or Boltzmann's constant [1.38 x 1(T 16 

g • cm 2 /(s 2 °K)], ML 2 lt 2 T or EIT y or 

mass transfer coefficient, Lit 
kf forward reaction rate constant, (mol/ 

£ 3 ) ! " w /r, n = order of reaction, various 

units 

ki mass transfer coefficient for species i, 
Lit 

k r reverse reaction rate constant, (mol/ 
L 3 ) 1_fl /r, n = order of reaction, various 
units 

K distribution coefficient, or hydrodynam- 
ic constant, dimensionless, or K—Cnbl 
ti y convenient dimensionless group in 
the dual-mode sorption model 

K c cake permeability, L 2 





uinusion consiani reiaung ainusiviiy to 


Q 1 


excess p 




concentration ior species z, l, /moi x or 


Irlt 






Q 


dimensio 




equilibrium constant (for extraction 


Qcr 


critical i 




reaction) , dimensionless or various units 




formatioi 




distribution coefficient for species i, 


Gv 


volumeti 




aja\ y dimensionless 


r 


radius, c 


Kreq 


equilibrium constant for stripping or re- 


r d 


drum rac 




extraction reaction, dimensionless or 


n 


inside ra 




various units 


To 


outside 


5 ■ 


solubility product, various units, or 




radius, I 




sorption coefficient, mol/Zr or MIL* 


r P 


pore rad 


K, 


tortuosity correction factor for bubble 


r, 


tube rad: 




point, dimensionless 


R 


gas cons 


i 
I 


membrane thickness, L 




of solutt 


L 


membrane length, or height of a column 




less 




extractor, L 


R c 


cake re: 


LRV 


filter log reduction value, dimensionless 




Lr l 


m 


molal concentration, mol/Af 


R c 


specific 


m, 


distribution coefficient for species i, 




ness, Lr 




(lila), dimensionless 


R'c 


specific 




mass, m 




M 


JV1 


molarity, rnoi/jL vg~moie/iuer or Kg- 


Re 


electrica 




mole/m ) 


Re 


Reynold 


M 


molecular weight, A//mol 


R f 


reaction 


n 


number of moles 




terfacial 


n p 


number of pores per unit membrane 




m 2 s) 




area, L z 


Ri 


rejectior 


N 


* /"lit /*! 1 /* 

number of hollow fibers, or number of 




tion par 




emulsion globules 


Rm 


membra 


XT 

N 


normal (g-equivalent/liter or kg-equiva- 




T -1 

or L 




lent/m ) 


Rmf 


final me 


AT 

N 


Avogadro s number, 0.023 x io J ■ 


Rmi 


initial ir 




molecules/mol 


Rr 


reaction 


A/7 


pressure drop, MILt 2 ox p 




area, m» 


p 


pressure or partial pressure, MILt 2 or p 


S 


cake coi 


p 


permeability coefficient, various units, 




less 




or fractional penetration, dimensionless 


S 


entropy . 


Pe 


Peclet number, vdlD, or v e LID u , di- 




for flu: 




mensionless 




geometi 


q 


heat flux, M/r or EIL z t 


S c 


specific 


Q 


flow rate or permeation rate, mol/r or 




solids v 




M/f, or volumetric flow rate, L 3 lt 


Sc 


Schmidi 




volumetric flow rate of the external 


Sh 

Oil 


Sherwo* 




phase, L 3 lt 




less 


Qi 


volumetric flow rate of the internal 


s„ 


specific 




phase, l?lt 




solids v 


Qm 


volumetric flow rate of the membrane 


t 


time, / 




phase, L?lt 


ha 


half-life 



General Notation xix 



g diffusivity to 
/, L 5 /mol t or 

for extraction 
>r various units 
for species /, 

stripping or re- 
lensionless or 

ous units, or 
II? or MIL 3 
tor for bubble 



;ht of a column 

dimensionless 

JM 

for species i, 



le/liter or kg- 
I 

nit membrane 
or number of 
or kg-equiva- 

023 x 10 23 



re, M/Lt 2 or p 
various units, 
dimensionless 
■ v e LID Uy di- 



rate, mol/f or 
rate, l?lt , 
the external 

: the internal 

lie membrane 



Q 

Qcr 
Cv 

r 
r d 

n 



r, 
R 



R c 
R, 



Re 



Ri 

Rr 
S 

S 



Sc 
Sh 



t 

ha 



excess particle flux in boundary layer, 
L 2 lt 

dimensionless excess particle flux 
critical excess particle flux for cake 
formation, dimensionless 
volumetric flow rate, L?it 
radius, or radial coordinate, L 
drum radius , L 
inside radius, L 

outside radius, or emulsion globule 
radius, L 
pore radius, L 
tube radius, L 

gas constant, ML 2 lt 2 T mol, or rejection 
of solute, or recycle ratio, dimension- 
less 

cake resistance, MIL 2 t (or ptIL), or 
L- 1 

specific cake resistance per unit thick- 
ness, Lr 2 

specific cake resistance (mass basis), LI . 
M 

electrical resistance, VIA 
Reynolds number, dvp/r), dimensionless 
reaction rate for extraction per in- 
terfacial area, mol f L 2 t (or kg-equivalent/ 
m 2 s) 

rejection of solute species /, or intercep- 
tion parameter, dimensionless 
membrane resistance, MIL 2 t (or ptIL), 
orZT 1 

final membrane resistance, L" 1 
initial membrane resistance, L" 1 
reaction rate for stripping per interfacial 
area, mol/L 2 f (or kg-equivalent/m 2 s) 
cake compressibility factor, dimension- 
less 

entropy, ML 2 lt 2 TojEIT, or shape factor 
for flux correction in a cylindrical 
geometry, dimensionless 
specific surface area (pore surface area/ 
solids volume) in cake, Lr 1 
Schmidt number, 17/pD, dimensionless 
Sherwood number, fc/L/D, dimension- 
less 

specific surface area (pore surface area/ 
solids volume) in membrane, L~ l 
time, t 
half-life, t 



T 

T b 

T c 

V 

v, 
V 



Vs 

Vi 
w 

w x 
W 

W e 

We 



yi 



transport number of ionic species i, di- 
mensionless 
temperature, T 
normal boiling point, T 
critical temperature, T 
glass transition temperature, T 
absolute mobility of species 1, t moVM 
or L 2 moVEt 
velocity, Lit 

velocity of the external phase, Lit 
inertial lift velocity for constricted tube 
or channel, Lit 

inertial lift velocity for unconstricted 
tube or channel, Lit 
superficial face velocity, Lit 
volume, l? 

volume of the external phase, L 3 

fractional free volume, dimensionless 

total volume of the internal phase, L 3 

volume of the membrane phase, L 3 

suspension volume, I? 

partial molar volume of species z, l?l 

mol 

cake mass per unit membrane area, Ml 
L 2 

mass fraction of species 1, dimension- 
less 

penetrant mass fraction, dimensionless 
channel width, or drum width, or width 
of membrane sheet, L 
total consumption amount of the reagent 
in the external phase, kg-equivalent 
Weber number, w 2 d 3 p e ly y dimension- 
less 

total consumption amount of the reagent 
in the internal phase, kg-equivalent 
critical distance for cake formation, L 
mole fraction of species /, dimension- 
less 

mole fraction of species / in permeate, 
dimensionless 

valence, dimensionless, or z coordinate, 
L 

valence of species 1, dimensionless 



Greek Letters 



selectivity, or mobility ratio (D AB c T l 
D a Cao) for facilitated transport, or 
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flow distribution parameter, dimension- 
less 

separation factor, (y/'/yJVW^Oi = 



y 

r 

r 



r, 



5c 



7?0 
Vd 

Vi 

Vm 
Vs 

e 



IT 
P 

Po 

Pe 
Pm 
Ps 

a 



c h Ph Wj, or,-, etc., dimensionless 
difference in polymer thermal expansion 
coefficients above and below glass tran- A 
sition, T" 1 (or °K _I ) /i 
matrix model coefficient in diffusion ex- 
pression, I? polymer/L 3 (STP) [or cm 3 v 
polymer/cm 3 (STP)] v, 
difference in polymer compressibility 
above and below glass transition, Lt 2 /M f 
or p~ x (or atm" 1 ) 

surface or interfacial tension, M/t 2 or 
Ell?, or plasticization parameter in the 
free-volume model, dimensionless 
activity coefficient of species i, dimen- 
sionless 
shear rate, t~ l 

shear rate at edge of boundary layer, t~ l 
tubesheet length, L 
inlet ratio of the internal reagent equiv- 
alents to the feed solute equivalents, di- a c 
mensionless 

ratio of the internal reagent equivalents 
to the external phase solute equivalents cr 0 
for stage i, dimensionless 
polarization layer thickness, or stagnant 
film thickness, L, or solubility param- 
eter (with subscript), (E/L 3 ) 1 ' 2 r 
cake thickness, L 

porosity, or void fraction, or in- 
homogeneity factor, dimensionless r c 
suspension porosity, dimensionless r m 
inverse Damkohler number, Da -1 , di- r w 
mensionless 

viscosity, M/Lt, or efficiency, dimen- Tho 
sionless 

pure fluid viscosity, M/Lt <f> 
single fiber efficiency due to diffusion, 
dimensionless <fi s 
viscosity of the external phase, M/Lt 
single fiber efficiency due to intercep- <f> w 
tion, dimensionless 

viscosity of the membrane phase, M/Lt i// 
solvent viscosity, M/Lt, or single col- <fc 
lector efficiency, dimensionless 
diffusion time lag, or time constant, /, <o 
or stage cut, or angle, dimensionless ft 



angle subtended by submerged drum, 
dimensionless 

ratio of species diameter to pore dia- 
meter, or parameter of the advancing 
front model, dimensionless 
reduced filter coefficient, dimensionless 
chemical potential, ML 2 It 2 mol or El 
mol 

kinematic viscosity, L 2 lt 
number of ions per molecule of electro- 
lyte i 

current utilization, or flow redistribu- 
tion decrement parameter, dimension- 
less 

osmotic pressure, M/Lt 2 or p 

density or mass concentration, MIL? 

pure fluid density, Mil? 

density of the external phase, MIL? 

density of the membrane phase, Mil? 

cake solids density, Mil? 

reflection coefficient, or inhomogeneity 

factor, dimensionless 

matrix model coefficient in solubility 

expression, I? polymer/L 3 (STP) [cm 3 

polymer/cm 3 (STP)] 

infinite dilution solubility coefficient in 

matrix model, I? (STP)/(L 3 polymer • 

M/Lt 2 ) (or cm 3 (STP)/(cm 3 polymer • 

atm)] 

tortuosity factor, or number of pore 
volumes of fluid passed through filter, 
dimensionless 

cake fouling time constant, t 
membrane fouling time constant, / 
wall shear stress for constricted tube or 
channel, M/Lt 2 

wall shear stress for unconstricted tube 
or channel, M/Lt 2 

volume fraction, or solids volume frac- 
tion, dimensionless 

suspension solids volume fraction, di- 
mensionless 

solids volume fraction at cake surface, 
dimensionless 
electric potential, V 
solvent association factor, dimension- 
less 

stirring rate or impeller speed, rpm 
angular velocity, t" ] 



Diacritical Mai 

per mol< 
per unit 
average 
time rat 



Superscripts 

I value ir 
stream s 
of the n 
external 

II value in 
value oi 
low-pre: 
value in 
sence o: 

0 standarc 
* ideal ca 



nerged drum, 

to pore dia- 
hc advancing 

ss 

dimensionless 
t 2 mol or £/ 



Diacritical Marks 



per mole 
per unit mass 
average value 
time rate of change 



ule of electro- 



I 



>w redistribu- \ Superscripts 
r, dimension- \ 



or p 

ition, Mil? 

lase, Mil? 
phase, MIL 1 

nhomogeneity 

: in solubility 
. 3 (ST?) [cm 3 

coefficient in 
(Z. 3 polymer ■ 
m 3 polymer • 

nber of pore 
. through Filter, 



it, / 

onstant, t 
tricted tube or 



Districted tube 
s volume frac- 
: fraction, di- 
cake surface, 



value in the feed stream, on the up- 
stream side or on the high-pressure side 
of the membrane, or value in the phase 
external to the membrane 
value in permeate, product or extract, or 
value on the downstream side or on the 
low-pressure side of the membrane 
value in the membrane indicated by ab- 
sence of superscript 
standard reference state 
idea] case 



Subscripts 

0 initial value 

A, B particular species 

b • bulk 

c cake 

/ feed 

r general species index or solute species i 

in inlet 

int interface 

J , species j 

/ liquid 

Im logarithmic mean 

tn membrane 

out outlet 

p product, permeate, permeant, or poly- 
mer 

r retentate or reject 

s solution or shell side 

t tube or tube side 

T total 

v vapor 

w water or solvent 

x t y t z three coordinate directions 



•r, dimension- 



;peed, rpm 



